HEAD ======= >>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4 ======= >>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
Data
library(survival)# the package for durvival functions.
library(ggfortify) # Usefull for intepreting Survival functions plots.
## Loading required package: ggplot2
library(KMsurv)
Cdata<-read.table("C_data.txt",header = FALSE)
head(Cdata)
=======
library(survival)# the package for durvival functions.
## Warning: package 'survival' was built under R version 4.0.3
library(ggfortify) # Usefull for intepreting Survival functions plots.
## Warning: package 'ggfortify' was built under R version 4.0.3
## Loading required package: ggplot2
Cdata<-read.table("C_data.txt",header = FALSE) #loading our data
head(Cdata)
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
## V1 V2 V3 V4 V5 V6 V7 V8
## 1 204 1 4 0 1 2 66 50
## 2 374 1 5 0 0 2 66 43
## 3 184 1 13 1 1 2 57 43
## 4 60 1 16 1 0 2 67 40
## 5 449 1 21 1 1 2 62 25
## 6 224 0 22 0 0 0 75 94
<<<<<<< HEAD
colnames(Cdata) <-c("Patient_number", "Death_Censoring", "Time", "Treatment", "Sex", "Ascites", "Age", "Prothrombin")
=======
colnames(Cdata) <-c("Patient_number", "Death_Censoring", "Time", "Treatment", "Sex", "Ascites", "Age", "Prothrombin") #adding names
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
head(Cdata)
## Patient_number Death_Censoring Time Treatment Sex Ascites Age Prothrombin
## 1 204 1 4 0 1 2 66 50
## 2 374 1 5 0 0 2 66 43
## 3 184 1 13 1 1 2 57 43
## 4 60 1 16 1 0 2 67 40
## 5 449 1 21 1 1 2 62 25
## 6 224 0 22 0 0 0 75 94
<<<<<<< HEAD
Opgave 1
Surv(Cdata$Time, Cdata$Death_Censoring)[1:10]
## [1] 4 5 13 16 21 22+ 23 24+ 24+ 25+
sur.fit1 <- survfit(Surv(Time, Death_Censoring) ~ Treatment, data=Cdata)
=======
sur.fit1 <- survfit(Surv(Time, Death_Censoring) ~ Treatment, data=Cdata)# constructing Survival function
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
names(sur.fit1)
## [1] "n" "time" "n.risk" "n.event" "n.censor" "surv"
## [7] "std.err" "cumhaz" "std.chaz" "strata" "type" "logse"
## [13] "conf.int" "conf.type" "lower" "upper" "call"
<<<<<<< HEAD
plot(sur.fit1,xlab="Days",ylab="Overall survival probability")
autoplot(sur.fit1,xlab="Days",ylab="Overall survival probability")
#plot(sur.fit1,xlab="Days",ylab="Overall survival probability") #a normal plot
autoplot(sur.fit1,xlab="Days",ylab="Overall survival probability") # a clear plot
# Log-rank test
Log_rank_test<-survdiff(Surv(Time, Death_Censoring) ~Treatment , data = Cdata)
Log_rank_test # results
## Call:
## survdiff(formula = Surv(Time, Death_Censoring) ~ Treatment, data = Cdata)
##
## N Observed Expected (O-E)^2/E (O-E)^2/V
## Treatment=0 251 142 149 0.355 0.728
## Treatment=1 237 150 143 0.371 0.728
##
## Chisq= 0.7 on 1 degrees of freedom, p= 0.4
<<<<<<< HEAD
1 - pchisq(Log_rank_test$chisq, length(Log_rank_test$n) - 1)#p-value
## [1] 0.3936019
Opgave 2
=======1 - pchisq(Log_rank_test$chisq, length(Log_rank_test$n) - 1)#True p-value
## [1] 0.3936019
The log-rank test show no significans in Treatment method
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
#Looking at Treatment
Cox.fit1 <- coxph(Sur1~ strata(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
summary(Cox.fit1)
## Call:
## coxph(formula = Sur1 ~ strata(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Sex)1 0.517727 1.678209 0.127220 4.070 4.71e-05 ***
## factor(Ascites)1 0.439503 1.551935 0.180240 2.438 0.0148 *
## factor(Ascites)2 0.966760 2.629412 0.182675 5.292 1.21e-07 ***
## Age 0.048929 1.050145 0.006810 7.185 6.74e-13 ***
## Prothrombin -0.012629 0.987450 0.002916 -4.331 1.49e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1 1.6782 0.5959 1.3078 2.1535
## factor(Ascites)1 1.5519 0.6444 1.0901 2.2095
## factor(Ascites)2 2.6294 0.3803 1.8381 3.7614
## Age 1.0501 0.9522 1.0362 1.0643
## Prothrombin 0.9875 1.0127 0.9818 0.9931
##
## Concordance= 0.692 (se = 0.018 )
## Likelihood ratio test= 126.1 on 5 df, p=<2e-16
## Wald test = 128.9 on 5 df, p=<2e-16
## Score (logrank) test = 138.4 on 5 df, p=<2e-16
source("andersen.R")
test1<-survfit(Cox.fit1)
test1
## Call: survfit(formula = Cox.fit1)
##
## n events median 0.95LCL 0.95UCL
## Treatment=0 251 142 1866 1374 2376
## Treatment=1 237 150 1588 1185 1978
andersen.plot(test1) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD
## [1] 0.002803319 0.005646071 0.005646071 0.005646071 0.005646071 0.005646071
## [7] 0.008523343 0.008523343 0.008523343 0.011437862 0.017335724 0.017335724
## [13] 0.020345066 0.023395335 0.023395335 0.023395335 0.026499061 0.029643468
## [19] 0.029643468 0.029643468 0.036035501 0.036035501 0.039307517 0.039307517
## [25] 0.039307517 0.039307517 0.042649418 0.042649418 0.046063196 0.046063196
## [31] 0.046063196 0.046063196 0.049541013 0.049541013 0.049541013 0.053049853
## [37] 0.056566027 0.056566027 0.060164114 0.063845205 0.063845205 0.063845205
## [43] 0.067587279 0.071356145 0.071356145 0.071356145 0.071356145 0.075181874
## [49] 0.079202607 0.083279132 0.087385608 0.087385608 0.091568985 0.095774186
## [55] 0.095774186 0.095774186 0.095774186 0.099991040 0.099991040 0.104263418
## [61] 0.112937833 0.117394860 0.117394860 0.121898044 0.121898044 0.121898044
## [67] 0.121898044 0.121898044 0.121898044 0.126483786 0.126483786 0.126483786
## [73] 0.126483786 0.131091385 0.131091385 0.131091385 0.140583177 0.140583177
## [79] 0.145384081 0.150193795 0.150193795 0.150193795 0.155024697 0.155024697
## [85] 0.155024697 0.164843078 0.169809657 0.174879989 0.174879989 0.174879989
## [91] 0.179983715 0.185148874 0.185148874 0.185148874 0.190343605 0.190343605
## [97] 0.195754728 0.201215599 0.201215599 0.201215599 0.206791483 0.212387598
## [103] 0.212387598 0.212387598 0.212387598 0.218077148 0.218077148 0.223830989
## [109] 0.223830989 0.223830989 0.229677281 0.235643235 0.235643235 0.241738991
## [115] 0.241738991 0.241738991 0.247896302 0.247896302 0.254199298 0.254199298
## [121] 0.254199298 0.254199298 0.260596686 0.260596686 0.267079498 0.267079498
## [127] 0.273660243 0.273660243 0.273660243 0.280259522 0.280259522 0.280259522
## [133] 0.286921714 0.286921714 0.286921714 0.293635839 0.293635839 0.300422372
## [139] 0.300422372 0.300422372 0.300422372 0.300422372 0.300422372 0.307299696
## [145] 0.314189319 0.314189319 0.314189319 0.314189319 0.314189319 0.314189319
## [151] 0.314189319 0.321146855 0.328152411 0.335176681 0.335176681 0.342270564
## [157] 0.342270564 0.342270564 0.342270564 0.342270564 0.349485048 0.349485048
## [163] 0.356733711 0.356733711 0.356733711 0.364048567 0.364048567 0.364048567
## [169] 0.364048567 0.364048567 0.371491815 0.371491815 0.371491815 0.379164770
## [175] 0.379164770 0.386915944 0.386915944 0.394805438 0.402733066 0.402733066
## [181] 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066
## [187] 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066
## [193] 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066 0.402733066
## [199] 0.402733066 0.411236396 0.411236396 0.411236396 0.411236396 0.411236396
## [205] 0.411236396 0.411236396 0.411236396 0.419795452 0.437088202 0.437088202
## [211] 0.437088202 0.437088202 0.437088202 0.437088202 0.437088202 0.437088202
## [217] 0.437088202 0.437088202 0.437088202 0.437088202 0.437088202 0.446533320
## [223] 0.446533320 0.456102394 0.465779829 0.465779829 0.465779829 0.465779829
## [229] 0.475600891 0.485506021 0.495492038 0.495492038 0.495492038 0.495492038
## [235] 0.495492038 0.495492038 0.505788223 0.505788223 0.516238757 0.526874067
## [241] 0.537591327 0.537591327 0.537591327 0.548523094 0.548523094 0.559879460
## [247] 0.571271313 0.571271313 0.571271313 0.582703773 0.582703773 0.594325938
## [253] 0.594325938 0.594325938 0.594325938 0.594325938 0.594325938 0.594325938
## [259] 0.594325938 0.606511480 0.606511480 0.618765886 0.631204366 0.631204366
## [265] 0.631204366 0.631204366 0.643858883 0.643858883 0.643858883 0.643858883
## [271] 0.656704522 0.669715610 0.669715610 0.682993370 0.696337571 0.696337571
## [277] 0.696337571 0.696337571 0.696337571 0.696337571 0.710249295 0.710249295
## [283] 0.710249295 0.710249295 0.710249295 0.710249295 0.710249295 0.710249295
## [289] 0.724349453 0.724349453 0.724349453 0.724349453 0.724349453 0.738577046
## [295] 0.738577046 0.753446696 0.753446696 0.753446696 0.753446696 0.753446696
## [301] 0.753446696 0.753446696 0.753446696 0.768898259 0.768898259 0.768898259
## [307] 0.768898259 0.768898259 0.768898259 0.768898259 0.768898259 0.785192921
## [313] 0.785192921 0.785192921 0.785192921 0.801765910 0.818684550 0.818684550
## [319] 0.818684550 0.836023320 0.836023320 0.853639281 0.871369433 0.889301550
## [325] 0.889301550 0.889301550 0.907645633 0.907645633 0.926439496 0.926439496
## [331] 0.926439496 0.926439496 0.926439496 0.947256706 0.947256706 0.947256706
## [337] 0.947256706 0.947256706 0.947256706 0.947256706 0.947256706 0.947256706
## [343] 0.947256706 0.947256706 0.969892125 0.969892125 0.969892125 0.969892125
## [349] 0.969892125 0.969892125 0.969892125 0.969892125 0.969892125 0.994528565
## [355] 0.994528565 0.994528565 0.994528565 0.994528565 0.994528565 0.994528565
## [361] 1.020675491 1.020675491 1.020675491 1.020675491 1.020675491 1.049912169
## [367] 1.079375654 1.079375654 1.079375654 1.079375654 1.079375654 1.079375654
## [373] 1.079375654 1.079375654 1.079375654 1.146811358 1.146811358 1.181835055
## [379] 1.218516764 1.218516764 1.256951458 1.298382879 1.298382879 1.341388730
## [385] 1.341388730 1.341388730 1.386772376 1.433078818 1.480546087 1.480546087
## [391] 1.480546087 1.529820305 1.529820305 1.529820305 1.529820305 1.529820305
## [397] 1.529820305 1.529820305 1.529820305 1.529820305 1.529820305 1.529820305
## [403] 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801
## [409] 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801 1.588748801
## [415] 1.588748801 1.588748801 1.588748801 1.706570781 1.851428550 1.851428550
## [421] 1.851428550 1.851428550 2.082565012 2.082565012 2.082565012 2.082565012
## [427] 2.082565012 2.082565012 2.082565012 2.082565012 2.082565012 2.082565012
## [433] 0.000000000 0.000000000 0.003250690 0.006538576 0.009864885 0.009864885
## [439] 0.009864885 0.009864885 0.009864885 0.013277201 0.013277201 0.013277201
## [445] 0.013277201 0.013277201 0.020207038 0.020207038 0.020207038 0.020207038
## [451] 0.020207038 0.023702448 0.027228455 0.030806238 0.034435535 0.038074158
## [457] 0.038074158 0.038074158 0.038074158 0.041732456 0.041732456 0.045403818
## [463] 0.049098253 0.049098253 0.049098253 0.052808653 0.056530257 0.056530257
## [469] 0.056530257 0.056530257 0.056530257 0.056530257 0.060256972 0.060256972
## [475] 0.060256972 0.060256972 0.060256972 0.060256972 0.060256972 0.071604831
## [481] 0.071604831 0.071604831 0.071604831 0.071604831 0.071604831 0.071604831
## [487] 0.079371547 0.083414992 0.087462234 0.087462234 0.099766915 0.099766915
## [493] 0.099766915 0.099766915 0.099766915 0.099766915 0.103915026 0.108113806
## [499] 0.108113806 0.108113806 0.108113806 0.108113806 0.112339438 0.112339438
## [505] 0.116587624 0.116587624 0.120853383 0.125300922 0.134296069 0.138899626
## [511] 0.138899626 0.143534473 0.143534473 0.148228377 0.148228377 0.152996332
## [517] 0.152996332 0.152996332 0.152996332 0.152996332 0.152996332 0.157794374
## [523] 0.162621192 0.162621192 0.162621192 0.167501229 0.172394281 0.177313034
## [529] 0.177313034 0.177313034 0.182249823 0.187236669 0.187236669 0.187236669
## [535] 0.192250769 0.197293558 0.202362380 0.207474451 0.207474451 0.207474451
## [541] 0.207474451 0.207474451 0.207474451 0.207474451 0.212634116 0.212634116
## [547] 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116
## [553] 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116 0.212634116
## [559] 0.212634116 0.218014502 0.218014502 0.218014502 0.223448922 0.223448922
## [565] 0.223448922 0.223448922 0.234510155 0.234510155 0.234510155 0.234510155
## [571] 0.240106474 0.245734315 0.251402356 0.251402356 0.251402356 0.251402356
## [577] 0.251402356 0.251402356 0.262881729 0.268674343 0.274480410 0.280323394
## [583] 0.286199096 0.286199096 0.286199096 0.286199096 0.286199096 0.286199096
## [589] 0.292129315 0.298073929 0.304107082 0.304107082 0.304107082 0.310162759
## [595] 0.310162759 0.316318573 0.322505090 0.322505090 0.322505090 0.322505090
## [601] 0.322505090 0.322505090 0.322505090 0.329083016 0.329083016 0.329083016
## [607] 0.335775894 0.335775894 0.342588674 0.342588674 0.342588674 0.342588674
## [613] 0.349410916 0.349410916 0.356306234 0.356306234 0.363250084 0.370238570
## [619] 0.377324446 0.377324446 0.384440778 0.391586591 0.391586591 0.398786655
## [625] 0.406261565 0.413790884 0.421370627 0.428975324 0.428975324 0.436859628
## [631] 0.444804484 0.444804484 0.452830595 0.460901444 0.469010410 0.477281466
## [637] 0.485627288 0.494045343 0.494045343 0.494045343 0.494045343 0.502797105
## [643] 0.511776939 0.511776939 0.511776939 0.511776939 0.520903087 0.530109827
## [649] 0.530109827 0.539428699 0.548789215 0.558218793 0.558218793 0.558218793
## [655] 0.558218793 0.567719358 0.567719358 0.567719358 0.577406531 0.587144719
## [661] 0.587144719 0.587144719 0.587144719 0.587144719 0.587144719 0.587144719
## [667] 0.587144719 0.597190566 0.597190566 0.607392500 0.607392500 0.607392500
## [673] 0.607392500 0.607392500 0.618050787 0.618050787 0.618050787 0.618050787
## [679] 0.618050787 0.618050787 0.628915165 0.628915165 0.628915165 0.628915165
## [685] 0.628915165 0.628915165 0.640015400 0.651237007 0.662783323 0.674765510
## [691] 0.686841825 0.686841825 0.699071671 0.699071671 0.699071671 0.699071671
## [697] 0.711603645 0.724294283 0.724294283 0.737229093 0.750560072 0.764056687
## [703] 0.764056687 0.764056687 0.777823629 0.777823629 0.777823629 0.777823629
## [709] 0.777823629 0.777823629 0.792009479 0.806448636 0.806448636 0.806448636
## [715] 0.821296645 0.821296645 0.836526845 0.852054132 0.867775447 0.883773968
## [721] 0.883773968 0.899977247 0.916477011 0.933106732 0.950032324 0.950032324
## [727] 0.950032324 0.950032324 0.950032324 0.950032324 0.950032324 0.950032324
## [733] 0.950032324 0.950032324 0.950032324 0.950032324 0.950032324 0.968763551
## [739] 0.987794194 1.007279790 1.007279790 1.007279790 1.027323522 1.027323522
## [745] 1.027323522 1.027323522 1.048040964 1.048040964 1.048040964 1.048040964
## [751] 1.048040964 1.048040964 1.069921181 1.069921181 1.069921181 1.069921181
## [757] 1.069921181 1.069921181 1.069921181 1.069921181 1.069921181 1.094575178
## [763] 1.119831131 1.119831131 1.148143865 1.148143865 1.148143865 1.177604111
## [769] 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111
## [775] 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111 1.177604111
## [781] 1.209863317 1.243221810 1.243221810 1.243221810 1.243221810 1.243221810
## [787] 1.277550368 1.277550368 1.313206852 1.313206852 1.313206852 1.351565798
## [793] 1.351565798 1.351565798 1.351565798 1.351565798 1.351565798 1.351565798
## [799] 1.351565798 1.351565798 1.351565798 1.351565798 1.351565798 1.392826090
## [805] 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090
## [811] 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090
## [817] 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090 1.392826090
## [823] 1.446508714 1.446508714 1.446508714 1.446508714 1.446508714 1.446508714
## [829] 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976
## [835] 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976
## [841] 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976 1.520464976
## [847] 1.669782450 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895
## [853] 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895
## [859] 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895 1.840477895
res1<-residuals(Cox.fit1,type="deviance")
hist(res1)
<<<<<<< HEAD
boxplot(res1~Cdata$Death_Censoring)
boxplot(res1~Cdata$Death_Censoring)
#Looking at sex
Cox.fit2 <- coxph(Sur1~ factor(Treatment)+ strata(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
summary(Cox.fit2)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + strata(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.054861 1.056393 0.118279 0.464 0.6428
## factor(Ascites)1 0.455062 1.576270 0.179616 2.534 0.0113 *
## factor(Ascites)2 0.977054 2.656619 0.182538 5.353 8.67e-08 ***
## Age 0.048919 1.050135 0.006785 7.209 5.62e-13 ***
## Prothrombin -0.012850 0.987233 0.002937 -4.376 1.21e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.0564 0.9466 0.8378 1.3320
## factor(Ascites)1 1.5763 0.6344 1.1085 2.2414
## factor(Ascites)2 2.6566 0.3764 1.8576 3.7993
## Age 1.0501 0.9523 1.0363 1.0642
## Prothrombin 0.9872 1.0129 0.9816 0.9929
##
## Concordance= 0.703 (se = 0.018 )
## Likelihood ratio test= 125 on 5 df, p=<2e-16
## Wald test = 129.3 on 5 df, p=<2e-16
## Score (logrank) test = 140.6 on 5 df, p=<2e-16
source("andersen.R")
test2<-survfit(Cox.fit2)
test2
## Call: survfit(formula = Cox.fit2)
##
## n events median 0.95LCL 0.95UCL
## Sex=0 198 111 2455 1950 2994
## Sex=1 290 181 1373 1141 1729
andersen.plot(test2) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD
## [1] 0.000000000 0.003345221 0.003345221 0.006739532 0.006739532 0.006739532
## [7] 0.010207978 0.010207978 0.010207978 0.010207978 0.010207978 0.010207978
## [13] 0.013778983 0.013778983 0.017427731 0.017427731 0.017427731 0.017427731
## [19] 0.017427731 0.021084701 0.024789265 0.028503991 0.032298009 0.032298009
## [25] 0.032298009 0.032298009 0.032298009 0.036167088 0.036167088 0.036167088
## [31] 0.040056452 0.040056452 0.043970676 0.043970676 0.043970676 0.047938128
## [37] 0.047938128 0.047938128 0.051917872 0.051917872 0.056005352 0.056005352
## [43] 0.056005352 0.056005352 0.056005352 0.056005352 0.056005352 0.060109904
## [49] 0.060109904 0.064231724 0.064231724 0.064231724 0.068394312 0.072585660
## [55] 0.072585660 0.076812794 0.076812794 0.076812794 0.089582682 0.093887342
## [61] 0.098227289 0.098227289 0.098227289 0.102626326 0.102626326 0.102626326
## [67] 0.102626326 0.102626326 0.102626326 0.107108781 0.111612664 0.111612664
## [73] 0.111612664 0.116159004 0.116159004 0.116159004 0.130114809 0.130114809
## [79] 0.134930460 0.139757868 0.139757868 0.139757868 0.139757868 0.139757868
## [85] 0.139757868 0.139757868 0.139757868 0.139757868 0.139757868 0.139757868
## [91] 0.144647340 0.144647340 0.144647340 0.144647340 0.144647340 0.149551589
## [97] 0.154481057 0.154481057 0.154481057 0.159466378 0.159466378 0.159466378
## [103] 0.164489861 0.169554440 0.174655924 0.179818794 0.179818794 0.185011809
## [109] 0.185011809 0.185011809 0.185011809 0.185011809 0.185011809 0.185011809
## [115] 0.185011809 0.185011809 0.185011809 0.185011809 0.190304939 0.190304939
## [121] 0.190304939 0.190304939 0.190304939 0.190304939 0.190304939 0.190304939
## [127] 0.190304939 0.195726388 0.195726388 0.195726388 0.201225033 0.201225033
## [133] 0.201225033 0.201225033 0.206830590 0.206830590 0.206830590 0.206830590
## [139] 0.206830590 0.206830590 0.212498753 0.212498753 0.212498753 0.218251218
## [145] 0.224014976 0.224014976 0.229841369 0.235707790 0.241593267 0.241593267
## [151] 0.247531873 0.253548931 0.253548931 0.253548931 0.253548931 0.253548931
## [157] 0.253548931 0.253548931 0.259625289 0.259625289 0.265733659 0.271874296
## [163] 0.271874296 0.271874296 0.278161128 0.278161128 0.278161128 0.278161128
## [169] 0.278161128 0.278161128 0.278161128 0.278161128 0.278161128 0.278161128
## [175] 0.278161128 0.278161128 0.284764700 0.291380525 0.291380525 0.291380525
## [181] 0.291380525 0.291380525 0.291380525 0.291380525 0.291380525 0.291380525
## [187] 0.291380525 0.291380525 0.298197782 0.305053024 0.305053024 0.312181122
## [193] 0.312181122 0.319696199 0.327281618 0.334901911 0.334901911 0.334901911
## [199] 0.334901911 0.342926589 0.350976673 0.359089262 0.367255347 0.367255347
## [205] 0.367255347 0.375658814 0.375658814 0.384589807 0.384589807 0.384589807
## [211] 0.384589807 0.384589807 0.384589807 0.384589807 0.394043327 0.394043327
## [217] 0.394043327 0.394043327 0.403874432 0.403874432 0.403874432 0.403874432
## [223] 0.403874432 0.413933604 0.413933604 0.413933604 0.413933604 0.413933604
## [229] 0.413933604 0.413933604 0.424213502 0.424213502 0.424213502 0.424213502
## [235] 0.424213502 0.424213502 0.424213502 0.424213502 0.434650119 0.445333746
## [241] 0.445333746 0.445333746 0.456279726 0.456279726 0.456279726 0.456279726
## [247] 0.467431862 0.467431862 0.467431862 0.467431862 0.467431862 0.478635809
## [253] 0.478635809 0.478635809 0.478635809 0.478635809 0.490022711 0.502016289
## [259] 0.502016289 0.514141966 0.514141966 0.514141966 0.514141966 0.514141966
## [265] 0.514141966 0.514141966 0.514141966 0.514141966 0.514141966 0.514141966
## [271] 0.514141966 0.514141966 0.514141966 0.526358492 0.538649435 0.538649435
## [277] 0.538649435 0.538649435 0.538649435 0.551125907 0.551125907 0.551125907
## [283] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [289] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [295] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [301] 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090 0.563830090
## [307] 0.577806544 0.577806544 0.577806544 0.577806544 0.577806544 0.592307003
## [313] 0.592307003 0.592307003 0.607040920 0.607040920 0.607040920 0.607040920
## [319] 0.607040920 0.607040920 0.607040920 0.622298172 0.637667908 0.653237280
## [325] 0.653237280 0.653237280 0.653237280 0.653237280 0.669499364 0.669499364
## [331] 0.687825832 0.687825832 0.708460386 0.729552770 0.729552770 0.750997675
## [337] 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675
## [343] 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675 0.750997675
## [349] 0.776640900 0.776640900 0.776640900 0.776640900 0.776640900 0.776640900
## [355] 0.803716523 0.803716523 0.803716523 0.803716523 0.803716523 0.833312019
## [361] 0.833312019 0.833312019 0.833312019 0.833312019 0.833312019 0.865278809
## [367] 0.865278809 0.865278809 0.865278809 0.865278809 0.865278809 0.899320721
## [373] 0.899320721 0.899320721 0.899320721 0.937040991 0.937040991 0.937040991
## [379] 0.937040991 0.937040991 0.975922999 0.975922999 0.975922999 1.018994173
## [385] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [391] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [397] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [403] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [409] 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173 1.018994173
## [415] 1.118012656 1.118012656 1.118012656 1.234982876 1.366600814 1.366600814
## [421] 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814
## [427] 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814 1.366600814
## [433] 0.002757535 0.002757535 0.005545448 0.005545448 0.008356406 0.008356406
## [439] 0.008356406 0.008356406 0.008356406 0.014074454 0.019861505 0.019861505
## [445] 0.019861505 0.022812378 0.025800620 0.025800620 0.028806331 0.031842293
## [451] 0.031842293 0.031842293 0.038015818 0.038015818 0.041181305 0.044352707
## [457] 0.044352707 0.044352707 0.047557766 0.047557766 0.050815624 0.054110208
## [463] 0.054110208 0.054110208 0.054110208 0.057430200 0.060757638 0.060757638
## [469] 0.064088457 0.064088457 0.064088457 0.067489583 0.067489583 0.067489583
## [475] 0.070931692 0.074391656 0.074391656 0.074391656 0.074391656 0.085107091
## [481] 0.088790926 0.088790926 0.092511874 0.092511874 0.092511874 0.092511874
## [487] 0.100137950 0.100137950 0.104080641 0.108053817 0.108053817 0.108053817
## [493] 0.112066311 0.116157297 0.116157297 0.116157297 0.120291324 0.124467166
## [499] 0.124467166 0.124467166 0.124467166 0.124467166 0.124467166 0.124467166
## [505] 0.128676372 0.128676372 0.132899856 0.137320430 0.141775934 0.146266732
## [511] 0.146266732 0.150782295 0.150782295 0.155326791 0.159929717 0.164567937
## [517] 0.164567937 0.173917460 0.178633284 0.183422929 0.183422929 0.188244562
## [523] 0.193090448 0.197980529 0.197980529 0.202928139 0.212933611 0.212933611
## [529] 0.212933611 0.218071708 0.223290566 0.223290566 0.228556180 0.233835838
## [535] 0.233835838 0.233835838 0.233835838 0.239180556 0.239180556 0.239180556
## [541] 0.239180556 0.239180556 0.244594070 0.250089607 0.255672777 0.261327717
## [547] 0.261327717 0.261327717 0.267073765 0.267073765 0.267073765 0.267073765
## [553] 0.267073765 0.267073765 0.272992839 0.272992839 0.278969209 0.278969209
## [559] 0.285010810 0.285010810 0.285010810 0.291064642 0.291064642 0.291064642
## [565] 0.297160780 0.297160780 0.303291099 0.309467267 0.309467267 0.315691443
## [571] 0.321969831 0.328280967 0.328280967 0.328280967 0.328280967 0.328280967
## [577] 0.328280967 0.328280967 0.334634071 0.334634071 0.334634071 0.341030821
## [583] 0.341030821 0.341030821 0.347459936 0.353901273 0.353901273 0.360380309
## [589] 0.366939999 0.373514305 0.373514305 0.373514305 0.373514305 0.373514305
## [595] 0.380178409 0.386886460 0.386886460 0.393624672 0.393624672 0.393624672
## [601] 0.393624672 0.393624672 0.400609212 0.407782608 0.407782608 0.415046290
## [607] 0.422365126 0.429802940 0.429802940 0.429802940 0.437342053 0.437342053
## [613] 0.444947615 0.444947615 0.452630371 0.452630371 0.460363210 0.468142839
## [619] 0.476022164 0.476022164 0.476022164 0.476022164 0.476022164 0.476022164
## [625] 0.483994395 0.483994395 0.483994395 0.483994395 0.483994395 0.492018105
## [631] 0.500093597 0.500093597 0.500093597 0.500093597 0.500093597 0.508238655
## [637] 0.516444093 0.516444093 0.516444093 0.516444093 0.533070772 0.541536117
## [643] 0.550180925 0.550180925 0.550180925 0.550180925 0.550180925 0.558938400
## [649] 0.558938400 0.567780493 0.567780493 0.576653625 0.576653625 0.585578905
## [655] 0.585578905 0.594577246 0.603650668 0.603650668 0.612858741 0.622104754
## [661] 0.631485851 0.640926781 0.640926781 0.640926781 0.640926781 0.640926781
## [667] 0.640926781 0.650667580 0.660529455 0.670502698 0.670502698 0.670502698
## [673] 0.680756087 0.680756087 0.680756087 0.691163264 0.691163264 0.701898395
## [679] 0.701898395 0.701898395 0.712689313 0.723583018 0.723583018 0.723583018
## [685] 0.723583018 0.723583018 0.734995825 0.746515007 0.746515007 0.746515007
## [691] 0.758324270 0.758324270 0.770255405 0.782427450 0.794741332 0.794741332
## [697] 0.807220965 0.819830004 0.832639334 0.845602503 0.858900126 0.872334656
## [703] 0.885990471 0.899793250 0.913836031 0.913836031 0.913836031 0.913836031
## [709] 0.913836031 0.913836031 0.928596975 0.928596975 0.943585879 0.943585879
## [715] 0.943585879 0.943585879 0.958987437 0.974640083 0.990455404 1.006504742
## [721] 1.022724644 1.039076112 1.055676826 1.072385926 1.089344418 1.106801640
## [727] 1.106801640 1.125021243 1.125021243 1.125021243 1.125021243 1.125021243
## [733] 1.125021243 1.125021243 1.125021243 1.144058979 1.144058979 1.163731681
## [739] 1.163731681 1.183675937 1.183675937 1.183675937 1.204101833 1.204101833
## [745] 1.204101833 1.204101833 1.204101833 1.225177632 1.246693520 1.246693520
## [751] 1.246693520 1.269109977 1.291888682 1.291888682 1.291888682 1.291888682
## [757] 1.291888682 1.291888682 1.316739239 1.316739239 1.316739239 1.342240244
## [763] 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244
## [769] 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244 1.342240244
## [775] 1.342240244 1.342240244 1.369335341 1.369335341 1.369335341 1.369335341
## [781] 1.369335341 1.397623227 1.397623227 1.397623227 1.397623227 1.426808943
## [787] 1.426808943 1.426808943 1.457454480 1.457454480 1.457454480 1.457454480
## [793] 1.488957202 1.488957202 1.488957202 1.488957202 1.488957202 1.488957202
## [799] 1.523228454 1.523228454 1.523228454 1.523228454 1.523228454 1.523228454
## [805] 1.523228454 1.523228454 1.523228454 1.560899224 1.560899224 1.600963914
## [811] 1.642737790 1.642737790 1.642737790 1.687062381 1.687062381 1.687062381
## [817] 1.687062381 1.687062381 1.735525345 1.784811988 1.835120582 1.835120582
## [823] 1.887003831 1.940567935 1.940567935 1.940567935 1.940567935 1.940567935
## [829] 2.006758591 2.006758591 2.006758591 2.006758591 2.006758591 2.006758591
## [835] 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732
## [841] 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732 2.087886732
## [847] 2.087886732 2.247980143 2.247980143 2.247980143 2.247980143 2.247980143
## [853] 2.247980143 2.247980143 2.599737118 2.599737118 2.599737118 2.599737118
## [859] 2.599737118 2.599737118 2.599737118 2.599737118 2.599737118 2.599737118
res2<-residuals(Cox.fit2,type="deviance")
hist(res2)
<<<<<<< HEAD
boxplot(res2~Cdata$Death_Censoring)
boxplot(res2~Cdata$Death_Censoring)
# Looking at Ascites
Cox.fit3 <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ strata(Ascites)+ Age+ Prothrombin, data = Cdata)
summary(Cox.fit3)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + strata(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.092631 1.097057 0.118711 0.780 0.43521
## factor(Sex)1 0.481967 1.619256 0.126547 3.809 0.00014 ***
## Age 0.048120 1.049297 0.006795 7.082 1.42e-12 ***
## Prothrombin -0.012599 0.987480 0.002929 -4.301 1.70e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.0971 0.9115 0.8693 1.3844
## factor(Sex)1 1.6193 0.6176 1.2636 2.0751
## Age 1.0493 0.9530 1.0354 1.0634
## Prothrombin 0.9875 1.0127 0.9818 0.9932
##
## Concordance= 0.646 (se = 0.02 )
## Likelihood ratio test= 75.12 on 4 df, p=2e-15
## Wald test = 70.1 on 4 df, p=2e-14
## Score (logrank) test = 70.43 on 4 df, p=2e-14
source("andersen.R")
test3<-survfit(Cox.fit3)
test3
## Call: survfit(formula = Cox.fit3)
##
## n events median 0.95LCL 0.95UCL
## Ascites=0 386 211 1979 1618 2277
## Ascites=1 54 39 1194 802 2187
## Ascites=2 48 42 686 224 1529
andersen.plot(test3) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD
## [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.002422910 0.002422910
## [7] 0.002422910 0.002422910 0.007342705 0.007342705 0.009818343 0.012299473
## [13] 0.012299473 0.017311969 0.019844167 0.022381833 0.022381833 0.022381833
## [19] 0.024936882 0.024936882 0.027499712 0.030075870 0.030075870 0.032669946
## [25] 0.035270579 0.037874232 0.037874232 0.037874232 0.040499620 0.043138885
## [31] 0.043138885 0.045803991 0.045803991 0.045803991 0.045803991 0.045803991
## [37] 0.051234459 0.053959659 0.053959659 0.056709190 0.056709190 0.056709190
## [43] 0.059476974 0.062246981 0.067836269 0.067836269 0.067836269 0.070640906
## [49] 0.070640906 0.073473549 0.076333898 0.076333898 0.076333898 0.076333898
## [55] 0.079228569 0.082130814 0.082130814 0.085048413 0.087975882 0.087975882
## [61] 0.090949553 0.096957268 0.099985028 0.103028855 0.109139889 0.109139889
## [67] 0.109139889 0.112228053 0.112228053 0.118449962 0.121585833 0.121585833
## [73] 0.124740883 0.127910783 0.131086222 0.131086222 0.134295773 0.137512089
## [79] 0.140741665 0.140741665 0.140741665 0.143980878 0.147245615 0.150524681
## [85] 0.153811349 0.153811349 0.157132578 0.160467711 0.167198221 0.167198221
## [91] 0.170591826 0.170591826 0.170591826 0.170591826 0.174035043 0.174035043
## [97] 0.174035043 0.174035043 0.174035043 0.174035043 0.174035043 0.177557058
## [103] 0.177557058 0.181106991 0.181106991 0.184688635 0.188276315 0.188276315
## [109] 0.188276315 0.191893457 0.191893457 0.195545936 0.195545936 0.199220479
## [115] 0.202918422 0.202918422 0.206640238 0.210388559 0.214153663 0.217940133
## [121] 0.217940133 0.217940133 0.221759485 0.225583152 0.225583152 0.233296347
## [127] 0.237181790 0.241074706 0.244987967 0.248918660 0.252879175 0.256857068
## [133] 0.260841555 0.260841555 0.260841555 0.264850422 0.268867250 0.272931469
## [139] 0.272931469 0.277008479 0.281097733 0.285209776 0.289338169 0.293505342
## [145] 0.293505342 0.293505342 0.293505342 0.293505342 0.297817238 0.297817238
## [151] 0.302187933 0.306586088 0.311044871 0.315508708 0.319986281 0.319986281
## [157] 0.319986281 0.319986281 0.324531278 0.324531278 0.324531278 0.329130385
## [163] 0.333778961 0.333778961 0.338472789 0.343182554 0.347921280 0.352686126
## [169] 0.357475797 0.362277850 0.367093534 0.371936042 0.376814323 0.381700710
## [175] 0.386607441 0.386607441 0.391531514 0.396487204 0.396487204 0.401486385
## [181] 0.406493245 0.406493245 0.411533108 0.411533108 0.411533108 0.411533108
## [187] 0.416642840 0.421783672 0.421783672 0.427027876 0.432287963 0.437574903
## [193] 0.437574903 0.442910860 0.442910860 0.453714946 0.459186349 0.459186349
## [199] 0.464725859 0.464725859 0.470285325 0.475873942 0.481490529 0.481490529
## [205] 0.481490529 0.481490529 0.481490529 0.487209959 0.492988895 0.498820289
## [211] 0.498820289 0.504785833 0.504785833 0.510816417 0.510816417 0.516962041
## [217] 0.523119100 0.523119100 0.529304753 0.535538855 0.535538855 0.541833971
## [223] 0.541833971 0.541833971 0.548306163 0.548306163 0.554827461 0.561383167
## [229] 0.567992686 0.567992686 0.574625244 0.581315779 0.581315779 0.588074133
## [235] 0.594886993 0.601782986 0.601782986 0.608740305 0.615751400 0.622848256
## [241] 0.622848256 0.629999870 0.637258708 0.644540098 0.644540098 0.644540098
## [247] 0.644540098 0.652030508 0.652030508 0.659604587 0.659604587 0.667301465
## [253] 0.667301465 0.675122748 0.683036373 0.691009874 0.699068259 0.707188520
## [259] 0.715354708 0.723609767 0.731903802 0.740285562 0.748838785 0.748838785
## [265] 0.757639804 0.757639804 0.757639804 0.757639804 0.757639804 0.757639804
## [271] 0.757639804 0.757639804 0.757639804 0.757639804 0.766967351 0.776382646
## [277] 0.785932909 0.785932909 0.785932909 0.795662298 0.805630556 0.805630556
## [283] 0.805630556 0.815715614 0.825948817 0.836324102 0.836324102 0.836324102
## [289] 0.847059561 0.857909699 0.868808764 0.879793185 0.879793185 0.879793185
## [295] 0.891329087 0.891329087 0.903058135 0.903058135 0.914948511 0.926971810
## [301] 0.926971810 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300
## [307] 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300
## [313] 0.939203300 0.939203300 0.952862217 0.966757789 0.966757789 0.966757789
## [319] 0.966757789 0.981074755 0.995864031 0.995864031 1.010953971 1.010953971
## [325] 1.010953971 1.026685796 1.042530764 1.042530764 1.042530764 1.042530764
## [331] 1.042530764 1.059580899 1.076967045 1.076967045 1.076967045 1.076967045
## [337] 1.076967045 1.076967045 1.076967045 1.076967045 1.116142086 1.116142086
## [343] 1.136740519 1.136740519 1.136740519 1.158208569 1.180687663 1.180687663
## [349] 1.204215922 1.204215922 1.204215922 1.229160124 1.254409050 1.280031424
## [355] 1.280031424 1.306648637 1.333899572 1.333899572 1.333899572 1.333899572
## [361] 1.333899572 1.333899572 1.333899572 1.333899572 1.333899572 1.333899572
## [367] 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164
## [373] 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164 1.367097164
## [379] 1.430615609 1.498616602 1.498616602 1.571643693 1.655346574 1.655346574
## [385] 1.655346574 1.655346574 1.780452655 1.780452655 1.780452655 1.780452655
## [391] 1.780452655 1.780452655 1.780452655 1.780452655 1.780452655 0.000000000
## [397] 0.012819623 0.012819623 0.012819623 0.026034998 0.039566407 0.039566407
## [403] 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500
## [409] 0.053402500 0.068122762 0.068122762 0.068122762 0.068122762 0.068122762
## [415] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
## [421] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
## [427] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.100386010
## [433] 0.100386010 0.100386010 0.100386010 0.117547315 0.134850559 0.134850559
## [439] 0.134850559 0.134850559 0.153076521 0.171646666 0.171646666 0.171646666
## [445] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
## [451] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
## [457] 0.191372504 0.191372504 0.191372504 0.191372504 0.191372504 0.211774184
## [463] 0.211774184 0.211774184 0.211774184 0.211774184 0.211774184 0.211774184
## [469] 0.232795610 0.232795610 0.232795610 0.232795610 0.232795610 0.232795610
## [475] 0.254560972 0.276919329 0.276919329 0.276919329 0.276919329 0.276919329
## [481] 0.300661572 0.300661572 0.300661572 0.300661572 0.300661572 0.300661572
## [487] 0.300661572 0.300661572 0.324925481 0.324925481 0.324925481 0.324925481
## [493] 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784
## [499] 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784 0.378294064
## [505] 0.378294064 0.378294064 0.378294064 0.378294064 0.406697290 0.406697290
## [511] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
## [517] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
## [523] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
## [529] 0.406697290 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806
## [535] 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806
## [541] 0.436020806 0.436020806 0.466837198 0.466837198 0.466837198 0.466837198
## [547] 0.466837198 0.466837198 0.466837198 0.466837198 0.466837198 0.499883881
## [553] 0.499883881 0.499883881 0.499883881 0.534298268 0.534298268 0.534298268
## [559] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
## [565] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
## [571] 0.569580748 0.569580748 0.569580748 0.569580748 0.569580748 0.607417459
## [577] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
## [583] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
## [589] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
## [595] 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010
## [601] 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010 0.756226472
## [607] 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
## [613] 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
## [619] 0.756226472 0.756226472 0.818519693 0.818519693 0.818519693 0.818519693
## [625] 0.889304890 0.889304890 0.889304890 0.889304890 0.889304890 0.889304890
## [631] 0.889304890 0.968545149 0.968545149 0.968545149 0.968545149 1.060605655
## [637] 1.060605655 1.060605655 1.060605655 1.060605655 1.060605655 1.060605655
## [643] 1.060605655 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [649] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [655] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [661] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [667] 1.162895503 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [673] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [679] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [685] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [691] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [697] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [703] 1.288056766 1.288056766 1.288056766 1.472408461 1.472408461 1.472408461
## [709] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [715] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [721] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [727] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [733] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [739] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877
## [745] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877
## [751] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.986477110
## [757] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
## [763] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
## [769] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
## [775] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
## [781] 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110 1.986477110
## [787] 1.986477110 1.986477110 1.986477110 1.986477110 0.061692868 0.061692868
## [793] 0.061692868 0.061692868 0.061692868 0.074960750 0.074960750 0.088612338
## [799] 0.088612338 0.088612338 0.088612338 0.088612338 0.088612338 0.117433167
## [805] 0.132798388 0.132798388 0.132798388 0.132798388 0.148612854 0.148612854
## [811] 0.148612854 0.148612854 0.148612854 0.165103160 0.165103160 0.165103160
## [817] 0.181895734 0.199245046 0.199245046 0.199245046 0.199245046 0.199245046
## [823] 0.217183883 0.217183883 0.217183883 0.217183883 0.235394131 0.255492719
## [829] 0.276175004 0.276175004 0.276175004 0.297464973 0.297464973 0.297464973
## [835] 0.345150450 0.345150450 0.369758625 0.369758625 0.369758625 0.395949724
## [841] 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724
## [847] 0.395949724 0.395949724 0.395949724 0.422705551 0.422705551 0.451879735
## [853] 0.451879735 0.451879735 0.451879735 0.451879735 0.483048139 0.483048139
## [859] 0.483048139 0.483048139 0.515749199 0.515749199 0.515749199 0.515749199
## [865] 0.515749199 0.515749199 0.515749199 0.550363970 0.550363970 0.550363970
## [871] 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
## [877] 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
## [883] 0.550363970 0.550363970 0.630418634 0.630418634 0.630418634 0.630418634
## [889] 0.630418634 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [895] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [901] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [907] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [913] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [919] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [925] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [931] 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
## [937] 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
## [943] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924
## [949] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924
## [955] 0.775050924 0.775050924 0.837188229 0.837188229 0.837188229 0.837188229
## [961] 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229
## [967] 0.837188229 0.837188229 0.922115104 0.922115104 0.922115104 0.922115104
## [973] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [979] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [985] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [991] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [997] 0.922115104 0.922115104 0.922115104 0.922115104 1.025522965 1.025522965
## [1003] 1.025522965 1.025522965 1.137150675 1.137150675 1.137150675 1.137150675
## [1009] 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675
## [1015] 1.137150675 1.137150675 1.269696765 1.269696765 1.269696765 1.269696765
## [1021] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1027] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1033] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1039] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1045] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1051] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1057] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1063] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1069] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1075] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.434236318
## [1081] 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318
## [1087] 1.640072935 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1093] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1099] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1105] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1111] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1117] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1123] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.781716629
## [1129] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1135] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1141] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1147] 2.781716629
## [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [7] 0.000000000 0.002422910 0.002422910 0.002422910 0.002422910 0.007342705
## [13] 0.007342705 0.009818343 0.012299473 0.012299473 0.012299473 0.017311969
## [19] 0.017311969 0.019844167 0.022381833 0.022381833 0.022381833 0.024936882
## [25] 0.027499712 0.030075870 0.030075870 0.030075870 0.032669946 0.035270579
## [31] 0.037874232 0.037874232 0.037874232 0.037874232 0.040499620 0.043138885
## [37] 0.043138885 0.045803991 0.045803991 0.045803991 0.045803991 0.045803991
## [43] 0.051234459 0.051234459 0.053959659 0.053959659 0.053959659 0.056709190
## [49] 0.056709190 0.056709190 0.059476974 0.062246981 0.062246981 0.067836269
## [55] 0.067836269 0.070640906 0.070640906 0.070640906 0.073473549 0.076333898
## [61] 0.076333898 0.076333898 0.076333898 0.079228569 0.082130814 0.082130814
## [67] 0.085048413 0.087975882 0.087975882 0.090949553 0.096957268 0.099985028
## [73] 0.103028855 0.109139889 0.109139889 0.109139889 0.112228053 0.112228053
## [79] 0.118449962 0.118449962 0.121585833 0.121585833 0.124740883 0.127910783
## [85] 0.131086222 0.131086222 0.134295773 0.137512089 0.140741665 0.143980878
## [91] 0.147245615 0.150524681 0.153811349 0.157132578 0.160467711 0.167198221
## [97] 0.167198221 0.170591826 0.170591826 0.170591826 0.170591826 0.170591826
## [103] 0.174035043 0.174035043 0.174035043 0.174035043 0.174035043 0.174035043
## [109] 0.174035043 0.174035043 0.177557058 0.177557058 0.181106991 0.181106991
## [115] 0.184688635 0.188276315 0.188276315 0.191893457 0.191893457 0.195545936
## [121] 0.195545936 0.199220479 0.202918422 0.202918422 0.206640238 0.210388559
## [127] 0.214153663 0.217940133 0.217940133 0.217940133 0.221759485 0.225583152
## [133] 0.225583152 0.233296347 0.237181790 0.241074706 0.244987967 0.248918660
## [139] 0.252879175 0.256857068 0.260841555 0.260841555 0.264850422 0.268867250
## [145] 0.272931469 0.272931469 0.277008479 0.277008479 0.281097733 0.285209776
## [151] 0.289338169 0.293505342 0.293505342 0.293505342 0.293505342 0.293505342
## [157] 0.293505342 0.297817238 0.297817238 0.302187933 0.306586088 0.306586088
## [163] 0.311044871 0.315508708 0.319986281 0.319986281 0.319986281 0.324531278
## [169] 0.324531278 0.329130385 0.333778961 0.333778961 0.338472789 0.343182554
## [175] 0.343182554 0.343182554 0.347921280 0.352686126 0.357475797 0.362277850
## [181] 0.362277850 0.367093534 0.371936042 0.376814323 0.381700710 0.386607441
## [187] 0.391531514 0.396487204 0.396487204 0.396487204 0.401486385 0.406493245
## [193] 0.411533108 0.411533108 0.411533108 0.416642840 0.421783672 0.421783672
## [199] 0.427027876 0.432287963 0.437574903 0.437574903 0.442910860 0.442910860
## [205] 0.453714946 0.459186349 0.459186349 0.464725859 0.470285325 0.475873942
## [211] 0.481490529 0.481490529 0.481490529 0.481490529 0.487209959 0.492988895
## [217] 0.498820289 0.498820289 0.504785833 0.504785833 0.510816417 0.510816417
## [223] 0.510816417 0.516962041 0.523119100 0.523119100 0.529304753 0.535538855
## [229] 0.541833971 0.541833971 0.541833971 0.548306163 0.548306163 0.554827461
## [235] 0.561383167 0.567992686 0.574625244 0.581315779 0.581315779 0.588074133
## [241] 0.594886993 0.601782986 0.608740305 0.615751400 0.622848256 0.629999870
## [247] 0.637258708 0.644540098 0.644540098 0.644540098 0.644540098 0.652030508
## [253] 0.659604587 0.659604587 0.667301465 0.667301465 0.675122748 0.683036373
## [259] 0.691009874 0.699068259 0.707188520 0.715354708 0.723609767 0.731903802
## [265] 0.740285562 0.748838785 0.748838785 0.757639804 0.757639804 0.757639804
## [271] 0.757639804 0.757639804 0.757639804 0.757639804 0.757639804 0.766967351
## [277] 0.776382646 0.785932909 0.785932909 0.785932909 0.795662298 0.805630556
## [283] 0.805630556 0.805630556 0.815715614 0.825948817 0.836324102 0.836324102
## [289] 0.836324102 0.847059561 0.847059561 0.857909699 0.868808764 0.879793185
## [295] 0.879793185 0.879793185 0.891329087 0.891329087 0.891329087 0.903058135
## [301] 0.903058135 0.914948511 0.926971810 0.926971810 0.939203300 0.939203300
## [307] 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300 0.939203300
## [313] 0.939203300 0.939203300 0.939203300 0.952862217 0.966757789 0.966757789
## [319] 0.966757789 0.966757789 0.981074755 0.995864031 0.995864031 1.010953971
## [325] 1.010953971 1.010953971 1.026685796 1.042530764 1.042530764 1.042530764
## [331] 1.042530764 1.042530764 1.059580899 1.076967045 1.076967045 1.076967045
## [337] 1.076967045 1.076967045 1.076967045 1.076967045 1.076967045 1.076967045
## [343] 1.116142086 1.116142086 1.136740519 1.136740519 1.158208569 1.180687663
## [349] 1.180687663 1.204215922 1.204215922 1.204215922 1.229160124 1.254409050
## [355] 1.280031424 1.280031424 1.306648637 1.333899572 1.333899572 0.000000000
## [361] 0.000000000 0.000000000 0.000000000 0.000000000 0.012819623 0.012819623
## [367] 0.026034998 0.039566407 0.039566407 0.039566407 0.053402500 0.053402500
## [373] 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500 0.053402500
## [379] 0.068122762 0.068122762 0.068122762 0.068122762 0.068122762 0.083749980
## [385] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
## [391] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980
## [397] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.100386010
## [403] 0.100386010 0.100386010 0.100386010 0.100386010 0.100386010 0.117547315
## [409] 0.134850559 0.134850559 0.134850559 0.134850559 0.134850559 0.171646666
## [415] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
## [421] 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666 0.171646666
## [427] 0.171646666 0.171646666 0.171646666 0.191372504 0.191372504 0.191372504
## [433] 0.191372504 0.191372504 0.191372504 0.211774184 0.211774184 0.211774184
## [439] 0.211774184 0.211774184 0.211774184 0.211774184 0.232795610 0.232795610
## [445] 0.232795610 0.232795610 0.232795610 0.232795610 0.276919329 0.276919329
## [451] 0.276919329 0.276919329 0.300661572 0.300661572 0.300661572 0.300661572
## [457] 0.300661572 0.300661572 0.300661572 0.324925481 0.324925481 0.324925481
## [463] 0.324925481 0.324925481 0.350680784 0.350680784 0.350680784 0.350680784
## [469] 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784 0.350680784
## [475] 0.350680784 0.350680784 0.378294064 0.378294064 0.378294064 0.378294064
## [481] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
## [487] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
## [493] 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290 0.406697290
## [499] 0.406697290 0.406697290 0.406697290 0.436020806 0.436020806 0.436020806
## [505] 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806
## [511] 0.436020806 0.436020806 0.436020806 0.436020806 0.436020806 0.466837198
## [517] 0.466837198 0.466837198 0.466837198 0.466837198 0.466837198 0.466837198
## [523] 0.466837198 0.466837198 0.466837198 0.499883881 0.499883881 0.499883881
## [529] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
## [535] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268
## [541] 0.534298268 0.534298268 0.534298268 0.534298268 0.534298268 0.569580748
## [547] 0.569580748 0.569580748 0.569580748 0.569580748 0.607417459 0.647833203
## [553] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
## [559] 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203 0.647833203
## [565] 0.647833203 0.647833203 0.647833203 0.698345010 0.698345010 0.698345010
## [571] 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010 0.698345010
## [577] 0.698345010 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
## [583] 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472 0.756226472
## [589] 0.756226472 0.756226472 0.756226472 0.818519693 0.818519693 0.818519693
## [595] 0.818519693 0.889304890 0.889304890 0.889304890 0.889304890 0.889304890
## [601] 0.889304890 0.968545149 0.968545149 0.968545149 1.060605655 1.060605655
## [607] 1.060605655 1.060605655 1.060605655 1.060605655 1.060605655 1.162895503
## [613] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [619] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [625] 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503 1.162895503
## [631] 1.162895503 1.162895503 1.162895503 1.288056766 1.288056766 1.288056766
## [637] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [643] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [649] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [655] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [661] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [667] 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766 1.288056766
## [673] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [679] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [685] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [691] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [697] 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461 1.472408461
## [703] 1.472408461 1.472408461 1.700899877 1.700899877 1.700899877 1.700899877
## [709] 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877 1.700899877
## [715] 1.700899877 1.700899877 1.700899877 1.700899877 0.023953755 0.036281742
## [721] 0.048854856 0.061692868 0.061692868 0.061692868 0.061692868 0.061692868
## [727] 0.074960750 0.074960750 0.088612338 0.088612338 0.088612338 0.088612338
## [733] 0.088612338 0.088612338 0.102901922 0.117433167 0.132798388 0.132798388
## [739] 0.132798388 0.132798388 0.148612854 0.148612854 0.148612854 0.148612854
## [745] 0.148612854 0.165103160 0.165103160 0.165103160 0.165103160 0.181895734
## [751] 0.181895734 0.199245046 0.199245046 0.199245046 0.199245046 0.199245046
## [757] 0.217183883 0.217183883 0.217183883 0.217183883 0.235394131 0.255492719
## [763] 0.255492719 0.276175004 0.276175004 0.276175004 0.276175004 0.297464973
## [769] 0.297464973 0.297464973 0.320944085 0.345150450 0.369758625 0.369758625
## [775] 0.369758625 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724
## [781] 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724 0.395949724
## [787] 0.422705551 0.422705551 0.451879735 0.451879735 0.451879735 0.451879735
## [793] 0.451879735 0.483048139 0.483048139 0.483048139 0.483048139 0.515749199
## [799] 0.515749199 0.515749199 0.515749199 0.515749199 0.515749199 0.515749199
## [805] 0.515749199 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
## [811] 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970 0.550363970
## [817] 0.550363970 0.550363970 0.589208062 0.630418634 0.630418634 0.630418634
## [823] 0.630418634 0.630418634 0.673825328 0.673825328 0.673825328 0.673825328
## [829] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [835] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [841] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [847] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [853] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [859] 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328 0.673825328
## [865] 0.673825328 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
## [871] 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533 0.720528533
## [877] 0.720528533 0.720528533 0.720528533 0.775050924 0.775050924 0.775050924
## [883] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924
## [889] 0.775050924 0.775050924 0.775050924 0.775050924 0.775050924 0.837188229
## [895] 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229
## [901] 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229 0.837188229
## [907] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [913] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [919] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [925] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104
## [931] 0.922115104 0.922115104 0.922115104 0.922115104 0.922115104 1.025522965
## [937] 1.025522965 1.025522965 1.025522965 1.137150675 1.137150675 1.137150675
## [943] 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675 1.137150675
## [949] 1.137150675 1.137150675 1.269696765 1.269696765 1.269696765 1.269696765
## [955] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [961] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [967] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [973] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [979] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [985] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [991] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [997] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1003] 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765 1.269696765
## [1009] 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318 1.434236318
## [1015] 1.434236318 1.434236318 1.640072935 1.640072935 2.013928967 2.013928967
## [1021] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1027] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1033] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1039] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1045] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1051] 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967 2.013928967
## [1057] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1063] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1069] 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629 2.781716629
## [1075] 2.781716629 2.781716629 2.781716629
## [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.002422910
## [7] 0.002422910 0.002422910 0.002422910 0.002422910 0.009818343 0.012299473
## [13] 0.012299473 0.017311969 0.017311969 0.019844167 0.022381833 0.022381833
## [19] 0.024936882 0.030075870 0.037874232 0.037874232 0.045803991 0.051234459
## [25] 0.051234459 0.053959659 0.056709190 0.056709190 0.062246981 0.067836269
## [31] 0.067836269 0.067836269 0.070640906 0.087975882 0.096957268 0.109139889
## [37] 0.109139889 0.118449962 0.127910783 0.137512089 0.140741665 0.140741665
## [43] 0.153811349 0.170591826 0.170591826 0.170591826 0.174035043 0.174035043
## [49] 0.174035043 0.188276315 0.199220479 0.260841555 0.277008479 0.293505342
## [55] 0.293505342 0.306586088 0.319986281 0.324531278 0.343182554 0.343182554
## [61] 0.362277850 0.386607441 0.396487204 0.406493245 0.406493245 0.411533108
## [67] 0.464725859 0.481490529 0.481490529 0.498820289 0.498820289 0.510816417
## [73] 0.535538855 0.548306163 0.548306163 0.567992686 0.601782986 0.622848256
## [79] 0.652030508 0.757639804 0.757639804 0.847059561 0.891329087 0.903058135
## [85] 0.903058135 0.939203300 0.939203300 1.076967045 1.136740519 1.333899572
## [91] 0.000000000 0.000000000 0.000000000 0.000000000 0.012819623 0.026034998
## [97] 0.039566407 0.039566407 0.039566407 0.053402500 0.053402500 0.053402500
## [103] 0.053402500 0.053402500 0.053402500 0.068122762 0.068122762 0.068122762
## [109] 0.083749980 0.083749980 0.083749980 0.083749980 0.083749980 0.100386010
## [115] 0.100386010 0.100386010 0.117547315 0.134850559 0.134850559 0.134850559
## [121] 0.153076521 0.171646666 0.171646666 0.171646666 0.191372504 0.191372504
## [127] 0.211774184 0.211774184 0.232795610 0.232795610 0.254560972 0.276919329
## [133] 0.300661572 0.324925481 0.324925481 0.324925481 0.350680784 0.350680784
## [139] 0.350680784 0.378294064 0.406697290 0.436020806 0.436020806 0.436020806
## [145] 0.466837198 0.466837198 0.499883881 0.534298268 0.534298268 0.534298268
## [151] 0.534298268 0.569580748 0.569580748 0.607417459 0.647833203 0.647833203
## [157] 0.698345010 0.698345010 0.698345010 0.698345010 0.756226472 0.756226472
## [163] 0.756226472 0.818519693 0.818519693 0.889304890 0.968545149 1.060605655
## [169] 1.162895503 1.162895503 1.288056766 1.288056766 1.288056766 1.288056766
## [175] 1.288056766 1.472408461 1.472408461 1.472408461 1.700899877 1.700899877
## [181] 0.023953755 0.036281742 0.048854856 0.061692868 0.061692868 0.061692868
## [187] 0.074960750 0.074960750 0.088612338 0.088612338 0.088612338 0.088612338
## [193] 0.102901922 0.117433167 0.132798388 0.132798388 0.132798388 0.148612854
## [199] 0.148612854 0.165103160 0.181895734 0.199245046 0.217183883 0.235394131
## [205] 0.255492719 0.276175004 0.276175004 0.297464973 0.320944085 0.345150450
## [211] 0.345150450 0.369758625 0.395949724 0.422705551 0.451879735 0.483048139
## [217] 0.483048139 0.515749199 0.515749199 0.550363970 0.550363970 0.550363970
## [223] 0.550363970 0.550363970 0.589208062 0.630418634 0.630418634 0.673825328
## [229] 0.673825328 0.673825328 0.673825328 0.673825328 0.720528533 0.720528533
## [235] 0.720528533 0.775050924 0.775050924 0.775050924 0.775050924 0.837188229
## [241] 0.837188229 0.837188229 0.922115104 0.922115104 0.922115104 0.922115104
## [247] 0.922115104 0.922115104 0.922115104 1.025522965 1.025522965 1.137150675
## [253] 1.137150675 1.137150675 1.269696765 1.269696765 1.269696765 1.269696765
## [259] 1.269696765 1.269696765 1.269696765 1.434236318 1.640072935 2.013928967
## [265] 2.013928967 2.013928967 2.013928967 2.781716629 2.781716629 2.781716629
res3<-residuals(Cox.fit3,type="deviance")
hist(res3)
<<<<<<< HEAD
boxplot(res3~Cdata$Death_Censoring)
boxplot(res3~Cdata$Death_Censoring)
# Looking at Age
min(Cdata[,7])# youngest person is 17
## [1] 17
max(Cdata[,7])# oldest person is 80
## [1] 80
Age.groups=cut(Cdata$Age,breaks=c(-Inf,quantile(Cdata$Age ,prob=c(0.33,0.66)),Inf)) #devides the Age into groups
Cox.fit4 <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ strata(Age.groups)+ Prothrombin, data = Cdata)
summary(Cox.fit4)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## strata(Age.groups) + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.026418 1.026770 0.119345 0.221 0.8248
## factor(Sex)1 0.497157 1.644041 0.127598 3.896 9.77e-05 ***
## factor(Ascites)1 0.459189 1.582789 0.183471 2.503 0.0123 *
## factor(Ascites)2 1.046229 2.846896 0.184316 5.676 1.38e-08 ***
## Prothrombin -0.012052 0.988020 0.002876 -4.191 2.78e-05 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.027 0.9739 0.8126 1.2974
## factor(Sex)1 1.644 0.6083 1.2803 2.1112
## factor(Ascites)1 1.583 0.6318 1.1047 2.2677
## factor(Ascites)2 2.847 0.3513 1.9837 4.0857
## Prothrombin 0.988 1.0121 0.9825 0.9936
##
## Concordance= 0.641 (se = 0.019 )
## Likelihood ratio test= 73.49 on 5 df, p=2e-14
## Wald test = 78.88 on 5 df, p=1e-15
## Score (logrank) test = 84.54 on 5 df, p=<2e-16
source("andersen.R")
test4<-survfit(Cox.fit4)
test4
## Call: survfit(formula = Cox.fit4)
##
## n events median 0.95LCL 0.95UCL
## (-Inf,55.7] 161 68 3078 2467 3611
## (55.7,64] 169 106 1559 1263 2057
## (64, Inf] 158 118 852 572 1169
andersen.plot(test4) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD
## [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.005528288
## [7] 0.005528288 0.005528288 0.005528288 0.005528288 0.005528288 0.005528288
## [13] 0.005528288 0.011403634 0.017304927 0.017304927 0.017304927 0.017304927
## [19] 0.017304927 0.023338937 0.029409752 0.029409752 0.029409752 0.035588296
## [25] 0.035588296 0.035588296 0.041919831 0.041919831 0.048303041 0.048303041
## [31] 0.048303041 0.054918882 0.061553184 0.068375151 0.068375151 0.068375151
## [37] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
## [43] 0.068375151 0.075365321 0.075365321 0.075365321 0.082455877 0.082455877
## [49] 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877
## [55] 0.082455877 0.082455877 0.089674197 0.096968283 0.096968283 0.104306718
## [61] 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718 0.111695071
## [67] 0.111695071 0.111695071 0.119142818 0.126825522 0.126825522 0.126825522
## [73] 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522
## [79] 0.126825522 0.134984552 0.134984552 0.143227902 0.143227902 0.143227902
## [85] 0.151634521 0.151634521 0.151634521 0.151634521 0.151634521 0.151634521
## [91] 0.151634521 0.160282070 0.160282070 0.168962279 0.168962279 0.177720498
## [97] 0.177720498 0.177720498 0.186541816 0.195438275 0.195438275 0.195438275
## [103] 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [109] 0.195438275 0.204585450 0.214054867 0.214054867 0.214054867 0.214054867
## [115] 0.223603523 0.223603523 0.233381313 0.243333679 0.243333679 0.253381794
## [121] 0.253381794 0.263515034 0.263515034 0.263515034 0.273710658 0.273710658
## [127] 0.273710658 0.273710658 0.283960762 0.283960762 0.283960762 0.283960762
## [133] 0.283960762 0.283960762 0.283960762 0.294431829 0.294431829 0.305000280
## [139] 0.305000280 0.305000280 0.305000280 0.305000280 0.305000280 0.305000280
## [145] 0.305000280 0.315826690 0.315826690 0.315826690 0.326983499 0.326983499
## [151] 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499 0.338812149
## [157] 0.338812149 0.338812149 0.338812149 0.350978763 0.363224224 0.363224224
## [163] 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224
## [169] 0.363224224 0.363224224 0.363224224 0.376029567 0.376029567 0.376029567
## [175] 0.376029567 0.376029567 0.389106667 0.389106667 0.402380512 0.416010968
## [181] 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813
## [187] 0.429762813 0.444009105 0.444009105 0.444009105 0.458561726 0.473262284
## [193] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [199] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.489049109
## [205] 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109
## [211] 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109 0.505600045
## [217] 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045
## [223] 0.505600045 0.523987973 0.523987973 0.543062495 0.543062495 0.543062495
## [229] 0.543062495 0.543062495 0.543062495 0.543062495 0.563339269 0.563339269
## [235] 0.563339269 0.563339269 0.586028574 0.586028574 0.586028574 0.586028574
## [241] 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574
## [247] 0.586028574 0.586028574 0.586028574 0.612099303 0.612099303 0.612099303
## [253] 0.612099303 0.612099303 0.639620829 0.639620829 0.639620829 0.639620829
## [259] 0.670502874 0.670502874 0.670502874 0.702589472 0.702589472 0.737090896
## [265] 0.737090896 0.737090896 0.737090896 0.777173206 0.818773352 0.818773352
## [271] 0.818773352 0.861528837 0.906277419 0.906277419 0.906277419 0.906277419
## [277] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [283] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [289] 0.959570305 0.959570305 0.959570305 1.090085142 1.090085142 1.090085142
## [295] 1.090085142 1.090085142 1.090085142 1.090085142 1.090085142 1.090085142
## [301] 1.090085142 1.090085142 1.090085142 0.005350627 0.010822944 0.010822944
## [307] 0.010822944 0.016454369 0.016454369 0.016454369 0.022164713 0.022164713
## [313] 0.027889662 0.027889662 0.027889662 0.033763774 0.033763774 0.033763774
## [319] 0.033763774 0.033763774 0.033763774 0.039791244 0.039791244 0.039791244
## [325] 0.039791244 0.045954956 0.045954956 0.045954956 0.045954956 0.045954956
## [331] 0.052253850 0.052253850 0.052253850 0.058725586 0.058725586 0.058725586
## [337] 0.065226295 0.065226295 0.071800183 0.071800183 0.071800183 0.071800183
## [343] 0.078528422 0.078528422 0.085296520 0.085296520 0.085296520 0.092181968
## [349] 0.099100984 0.099100984 0.106076035 0.106076035 0.113136184 0.120399425
## [355] 0.120399425 0.127705587 0.127705587 0.135076428 0.142522708 0.142522708
## [361] 0.142522708 0.149993842 0.149993842 0.157531961 0.165150610 0.165150610
## [367] 0.172846201 0.180620008 0.188445918 0.188445918 0.196368910 0.196368910
## [373] 0.196368910 0.196368910 0.196368910 0.196368910 0.204695686 0.204695686
## [379] 0.204695686 0.204695686 0.204695686 0.213256293 0.213256293 0.213256293
## [385] 0.213256293 0.221921101 0.221921101 0.221921101 0.230668576 0.230668576
## [391] 0.239491574 0.248425299 0.248425299 0.248425299 0.248425299 0.248425299
## [397] 0.257512175 0.266684841 0.266684841 0.266684841 0.275887140 0.275887140
## [403] 0.275887140 0.275887140 0.285301563 0.294762418 0.304319540 0.304319540
## [409] 0.304319540 0.304319540 0.304319540 0.314131316 0.314131316 0.314131316
## [415] 0.324057279 0.334041522 0.334041522 0.334041522 0.344145476 0.344145476
## [421] 0.344145476 0.344145476 0.344145476 0.354489364 0.354489364 0.364955963
## [427] 0.375468287 0.375468287 0.386067530 0.396810622 0.396810622 0.396810622
## [433] 0.418872475 0.430217353 0.430217353 0.430217353 0.441771494 0.453445406
## [439] 0.453445406 0.465319253 0.465319253 0.465319253 0.477384731 0.477384731
## [445] 0.489634257 0.502028814 0.502028814 0.514632492 0.514632492 0.527340450
## [451] 0.540164985 0.540164985 0.540164985 0.540164985 0.540164985 0.553304597
## [457] 0.566657491 0.580311979 0.580311979 0.580311979 0.594305960 0.594305960
## [463] 0.594305960 0.594305960 0.594305960 0.609012152 0.623904225 0.623904225
## [469] 0.639143157 0.639143157 0.654807088 0.670718435 0.687301167 0.704071054
## [475] 0.704071054 0.721157987 0.721157987 0.738561528 0.756346779 0.756346779
## [481] 0.774488197 0.774488197 0.774488197 0.774488197 0.774488197 0.793145557
## [487] 0.812226345 0.812226345 0.812226345 0.831644405 0.831644405 0.851361126
## [493] 0.871448014 0.871448014 0.871448014 0.891810367 0.891810367 0.913268483
## [499] 0.913268483 0.913268483 0.913268483 0.913268483 0.913268483 0.913268483
## [505] 0.913268483 0.936798339 0.936798339 0.936798339 0.936798339 0.961841394
## [511] 0.987712859 0.987712859 0.987712859 1.013960311 1.040953046 1.040953046
## [517] 1.040953046 1.069941379 1.069941379 1.099543534 1.129531930 1.129531930
## [523] 1.160547040 1.192692166 1.192692166 1.228912605 1.228912605 1.228912605
## [529] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [535] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [541] 1.228912605 1.228912605 1.228912605 1.270816176 1.270816176 1.314046197
## [547] 1.314046197 1.359029219 1.404934747 1.404934747 1.404934747 1.404934747
## [553] 1.404934747 1.458863406 1.458863406 1.458863406 1.458863406 1.458863406
## [559] 1.458863406 1.458863406 1.458863406 1.520684154 1.520684154 1.590254427
## [565] 1.590254427 1.590254427 1.590254427 1.590254427 1.667260347 1.667260347
## [571] 1.667260347 1.667260347 1.747291122 1.747291122 1.747291122 1.747291122
## [577] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [583] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [589] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.926625374
## [595] 1.926625374 1.926625374 1.926625374 1.926625374 1.926625374 2.281783691
## [601] 2.281783691 2.281783691 2.281783691 2.281783691 2.281783691 2.281783691
## [607] 0.009633169 0.014551755 0.019546218 0.019546218 0.024599344 0.029678509
## [613] 0.029678509 0.045413686 0.045413686 0.045413686 0.050750288 0.050750288
## [619] 0.050750288 0.061572580 0.072742518 0.078415489 0.078415489 0.107873484
## [625] 0.107873484 0.107873484 0.107873484 0.126231259 0.126231259 0.138970374
## [631] 0.138970374 0.138970374 0.158792804 0.158792804 0.165618521 0.165618521
## [637] 0.172484282 0.186450350 0.193667471 0.200947081 0.230731258 0.230731258
## [643] 0.246241191 0.246241191 0.246241191 0.246241191 0.254091780 0.254091780
## [649] 0.269933185 0.269933185 0.295740183 0.295740183 0.295740183 0.304607658
## [655] 0.304607658 0.304607658 0.313541365 0.313541365 0.322590148 0.341155661
## [661] 0.341155661 0.350635123 0.350635123 0.350635123 0.360361864 0.360361864
## [667] 0.360361864 0.391154668 0.401655837 0.401655837 0.401655837 0.412274302
## [673] 0.412274302 0.433912349 0.433912349 0.433912349 0.433912349 0.433912349
## [679] 0.445124076 0.445124076 0.445124076 0.445124076 0.445124076 0.456660219
## [685] 0.456660219 0.456660219 0.468424656 0.468424656 0.480317909 0.480317909
## [691] 0.492387762 0.492387762 0.492387762 0.492387762 0.504593246 0.516909136
## [697] 0.516909136 0.516909136 0.529347066 0.529347066 0.529347066 0.580175913
## [703] 0.593186674 0.593186674 0.593186674 0.606344153 0.606344153 0.606344153
## [709] 0.619675763 0.619675763 0.647161102 0.647161102 0.661946165 0.677209373
## [715] 0.677209373 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700
## [721] 0.692791700 0.692791700 0.692791700 0.708979609 0.708979609 0.708979609
## [727] 0.742673869 0.760665334 0.760665334 0.760665334 0.779916835 0.839297067
## [733] 0.839297067 0.859999536 0.859999536 0.859999536 0.881788673 0.881788673
## [739] 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673
## [745] 0.881788673 0.881788673 0.881788673 0.905686340 0.905686340 0.905686340
## [751] 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340
## [757] 0.905686340 0.905686340 0.930593182 0.956024835 0.956024835 0.956024835
## [763] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [769] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [775] 0.982860461 1.011138349 1.011138349 1.011138349 1.072147856 1.072147856
## [781] 1.104113664 1.104113664 1.136775819 1.171345121 1.171345121 1.207040278
## [787] 1.207040278 1.207040278 1.207040278 1.246907275 1.246907275 1.287930715
## [793] 1.330411123 1.330411123 1.330411123 1.330411123 1.374604413 1.374604413
## [799] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943
## [805] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.521709021
## [811] 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021
## [817] 1.574890670 1.574890670 1.574890670 1.574890670 1.574890670 1.574890670
## [823] 1.574890670 1.685791093 1.685791093 1.750311642 1.750311642 1.750311642
## [829] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [835] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [841] 1.750311642 1.750311642 1.832400858 1.832400858 1.832400858 1.832400858
## [847] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [853] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [859] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [865] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 2.029441753
## [871] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [877] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [883] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.194031226
## [889] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226
## [895] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.530841239
## [901] 2.925483659 2.925483659 2.925483659
## [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [7] 0.000000000 0.005528288 0.005528288 0.005528288 0.005528288 0.005528288
## [13] 0.005528288 0.005528288 0.005528288 0.011403634 0.011403634 0.017304927
## [19] 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927
## [25] 0.017304927 0.023338937 0.029409752 0.029409752 0.029409752 0.029409752
## [31] 0.029409752 0.029409752 0.029409752 0.029409752 0.035588296 0.035588296
## [37] 0.041919831 0.041919831 0.048303041 0.048303041 0.048303041 0.048303041
## [43] 0.054918882 0.054918882 0.061553184 0.068375151 0.068375151 0.068375151
## [49] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
## [55] 0.068375151 0.068375151 0.075365321 0.075365321 0.082455877 0.082455877
## [61] 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877
## [67] 0.082455877 0.089674197 0.089674197 0.096968283 0.096968283 0.104306718
## [73] 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718 0.111695071
## [79] 0.111695071 0.111695071 0.111695071 0.111695071 0.111695071 0.119142818
## [85] 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522
## [91] 0.126825522 0.126825522 0.126825522 0.134984552 0.134984552 0.134984552
## [97] 0.143227902 0.143227902 0.143227902 0.143227902 0.151634521 0.151634521
## [103] 0.151634521 0.151634521 0.151634521 0.160282070 0.160282070 0.168962279
## [109] 0.168962279 0.168962279 0.168962279 0.168962279 0.177720498 0.177720498
## [115] 0.177720498 0.186541816 0.186541816 0.195438275 0.195438275 0.195438275
## [121] 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [127] 0.195438275 0.195438275 0.204585450 0.214054867 0.214054867 0.223603523
## [133] 0.223603523 0.223603523 0.233381313 0.233381313 0.243333679 0.253381794
## [139] 0.253381794 0.253381794 0.253381794 0.253381794 0.263515034 0.263515034
## [145] 0.263515034 0.273710658 0.273710658 0.273710658 0.273710658 0.273710658
## [151] 0.283960762 0.283960762 0.283960762 0.283960762 0.283960762 0.294431829
## [157] 0.305000280 0.305000280 0.305000280 0.305000280 0.315826690 0.326983499
## [163] 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499 0.338812149
## [169] 0.338812149 0.338812149 0.338812149 0.350978763 0.363224224 0.363224224
## [175] 0.363224224 0.363224224 0.363224224 0.376029567 0.376029567 0.376029567
## [181] 0.376029567 0.376029567 0.376029567 0.389106667 0.389106667 0.402380512
## [187] 0.402380512 0.416010968 0.429762813 0.429762813 0.429762813 0.429762813
## [193] 0.429762813 0.429762813 0.429762813 0.429762813 0.444009105 0.444009105
## [199] 0.458561726 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [205] 0.473262284 0.473262284 0.473262284 0.489049109 0.489049109 0.489049109
## [211] 0.489049109 0.489049109 0.489049109 0.505600045 0.505600045 0.505600045
## [217] 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045 0.523987973
## [223] 0.523987973 0.543062495 0.543062495 0.543062495 0.543062495 0.543062495
## [229] 0.543062495 0.563339269 0.563339269 0.563339269 0.563339269 0.586028574
## [235] 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574
## [241] 0.586028574 0.586028574 0.612099303 0.612099303 0.612099303 0.612099303
## [247] 0.639620829 0.639620829 0.639620829 0.670502874 0.702589472 0.702589472
## [253] 0.702589472 0.737090896 0.737090896 0.737090896 0.737090896 0.777173206
## [259] 0.818773352 0.861528837 0.906277419 0.906277419 0.906277419 0.959570305
## [265] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [271] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [277] 0.959570305 1.090085142 1.090085142 1.090085142 1.090085142 1.090085142
## [283] 1.090085142 1.090085142 0.000000000 0.005350627 0.010822944 0.010822944
## [289] 0.010822944 0.010822944 0.016454369 0.016454369 0.016454369 0.016454369
## [295] 0.016454369 0.022164713 0.022164713 0.027889662 0.027889662 0.033763774
## [301] 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774
## [307] 0.033763774 0.033763774 0.033763774 0.039791244 0.039791244 0.039791244
## [313] 0.039791244 0.039791244 0.039791244 0.045954956 0.045954956 0.045954956
## [319] 0.045954956 0.045954956 0.045954956 0.052253850 0.052253850 0.052253850
## [325] 0.058725586 0.058725586 0.058725586 0.058725586 0.058725586 0.065226295
## [331] 0.065226295 0.065226295 0.065226295 0.071800183 0.071800183 0.071800183
## [337] 0.078528422 0.078528422 0.085296520 0.085296520 0.085296520 0.092181968
## [343] 0.099100984 0.106076035 0.113136184 0.120399425 0.127705587 0.127705587
## [349] 0.127705587 0.127705587 0.142522708 0.142522708 0.142522708 0.142522708
## [355] 0.149993842 0.149993842 0.157531961 0.157531961 0.157531961 0.165150610
## [361] 0.180620008 0.188445918 0.188445918 0.188445918 0.188445918 0.188445918
## [367] 0.188445918 0.196368910 0.196368910 0.196368910 0.196368910 0.196368910
## [373] 0.204695686 0.204695686 0.204695686 0.204695686 0.213256293 0.213256293
## [379] 0.213256293 0.213256293 0.213256293 0.213256293 0.213256293 0.221921101
## [385] 0.221921101 0.230668576 0.239491574 0.248425299 0.248425299 0.248425299
## [391] 0.248425299 0.257512175 0.266684841 0.266684841 0.266684841 0.266684841
## [397] 0.266684841 0.266684841 0.266684841 0.275887140 0.275887140 0.275887140
## [403] 0.275887140 0.285301563 0.304319540 0.304319540 0.304319540 0.304319540
## [409] 0.304319540 0.304319540 0.304319540 0.314131316 0.314131316 0.314131316
## [415] 0.334041522 0.334041522 0.334041522 0.344145476 0.344145476 0.344145476
## [421] 0.344145476 0.344145476 0.344145476 0.344145476 0.344145476 0.354489364
## [427] 0.354489364 0.364955963 0.375468287 0.375468287 0.375468287 0.375468287
## [433] 0.375468287 0.396810622 0.396810622 0.418872475 0.430217353 0.430217353
## [439] 0.453445406 0.453445406 0.465319253 0.477384731 0.489634257 0.502028814
## [445] 0.514632492 0.540164985 0.540164985 0.540164985 0.540164985 0.540164985
## [451] 0.553304597 0.580311979 0.580311979 0.580311979 0.594305960 0.594305960
## [457] 0.594305960 0.594305960 0.623904225 0.639143157 0.670718435 0.704071054
## [463] 0.704071054 0.704071054 0.721157987 0.721157987 0.756346779 0.756346779
## [469] 0.756346779 0.774488197 0.774488197 0.774488197 0.774488197 0.774488197
## [475] 0.774488197 0.793145557 0.812226345 0.812226345 0.812226345 0.812226345
## [481] 0.831644405 0.871448014 0.871448014 0.871448014 0.871448014 0.913268483
## [487] 0.913268483 0.913268483 0.936798339 0.936798339 0.936798339 0.936798339
## [493] 0.936798339 0.961841394 0.987712859 0.987712859 0.987712859 1.069941379
## [499] 1.099543534 1.099543534 1.099543534 1.099543534 1.129531930 1.129531930
## [505] 1.192692166 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [511] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [517] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.270816176
## [523] 1.314046197 1.314046197 1.404934747 1.404934747 1.404934747 1.458863406
## [529] 1.458863406 1.458863406 1.458863406 1.458863406 1.458863406 1.520684154
## [535] 1.590254427 1.590254427 1.590254427 1.590254427 1.590254427 1.667260347
## [541] 1.667260347 1.667260347 1.667260347 1.747291122 1.747291122 1.747291122
## [547] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [553] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [559] 1.747291122 1.747291122 1.747291122 1.926625374 1.926625374 1.926625374
## [565] 1.926625374 1.926625374 1.926625374 2.281783691 0.009633169 0.014551755
## [571] 0.014551755 0.019546218 0.019546218 0.019546218 0.024599344 0.029678509
## [577] 0.029678509 0.034864290 0.040107303 0.045413686 0.045413686 0.050750288
## [583] 0.050750288 0.061572580 0.067117647 0.072742518 0.078415489 0.078415489
## [589] 0.084121378 0.089991681 0.095877420 0.101853535 0.107873484 0.107873484
## [595] 0.107873484 0.113924256 0.119989528 0.119989528 0.126231259 0.132564626
## [601] 0.132564626 0.138970374 0.138970374 0.138970374 0.158792804 0.165618521
## [607] 0.165618521 0.172484282 0.172484282 0.186450350 0.186450350 0.193667471
## [613] 0.193667471 0.200947081 0.208250442 0.223100168 0.230731258 0.238451311
## [619] 0.246241191 0.246241191 0.254091780 0.254091780 0.261998834 0.269933185
## [625] 0.269933185 0.295740183 0.295740183 0.304607658 0.313541365 0.313541365
## [631] 0.322590148 0.331761267 0.341155661 0.341155661 0.350635123 0.350635123
## [637] 0.350635123 0.350635123 0.360361864 0.360361864 0.370472540 0.380732678
## [643] 0.391154668 0.401655837 0.412274302 0.412274302 0.412274302 0.422970130
## [649] 0.422970130 0.422970130 0.433912349 0.433912349 0.433912349 0.433912349
## [655] 0.445124076 0.445124076 0.445124076 0.445124076 0.456660219 0.456660219
## [661] 0.456660219 0.456660219 0.468424656 0.468424656 0.468424656 0.480317909
## [667] 0.480317909 0.480317909 0.492387762 0.492387762 0.504593246 0.516909136
## [673] 0.516909136 0.516909136 0.529347066 0.529347066 0.541876644 0.554515396
## [679] 0.567272570 0.580175913 0.580175913 0.593186674 0.593186674 0.593186674
## [685] 0.606344153 0.606344153 0.606344153 0.619675763 0.633344103 0.647161102
## [691] 0.647161102 0.647161102 0.647161102 0.661946165 0.677209373 0.692791700
## [697] 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700
## [703] 0.692791700 0.708979609 0.708979609 0.708979609 0.725379668 0.725379668
## [709] 0.742673869 0.760665334 0.760665334 0.760665334 0.779916835 0.779916835
## [715] 0.799534319 0.819325249 0.839297067 0.859999536 0.859999536 0.881788673
## [721] 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673
## [727] 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340
## [733] 0.905686340 0.930593182 0.956024835 0.956024835 0.956024835 0.982860461
## [739] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [745] 1.011138349 1.011138349 1.040995345 1.072147856 1.072147856 1.104113664
## [751] 1.136775819 1.136775819 1.171345121 1.171345121 1.207040278 1.207040278
## [757] 1.207040278 1.207040278 1.207040278 1.246907275 1.246907275 1.287930715
## [763] 1.287930715 1.330411123 1.330411123 1.374604413 1.374604413 1.374604413
## [769] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.470281401
## [775] 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021 1.521709021
## [781] 1.574890670 1.574890670 1.629630555 1.685791093 1.685791093 1.685791093
## [787] 1.685791093 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [793] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [799] 1.750311642 1.750311642 1.832400858 1.832400858 1.832400858 1.832400858
## [805] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [811] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [817] 1.917615102 1.917615102 1.917615102 1.917615102 2.029441753 2.029441753
## [823] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [829] 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [835] 2.029441753 2.029441753 2.029441753 2.194031226 2.194031226 2.194031226
## [841] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226
## [847] 2.194031226 2.530841239 2.925483659 2.925483659 2.925483659 2.925483659
## [1] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [7] 0.000000000 0.000000000 0.005528288 0.005528288 0.005528288 0.005528288
## [13] 0.005528288 0.005528288 0.005528288 0.011403634 0.011403634 0.017304927
## [19] 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927 0.017304927
## [25] 0.017304927 0.017304927 0.017304927 0.029409752 0.029409752 0.029409752
## [31] 0.029409752 0.029409752 0.029409752 0.029409752 0.029409752 0.035588296
## [37] 0.035588296 0.035588296 0.041919831 0.041919831 0.041919831 0.048303041
## [43] 0.048303041 0.048303041 0.048303041 0.054918882 0.068375151 0.068375151
## [49] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
## [55] 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151 0.068375151
## [61] 0.068375151 0.075365321 0.075365321 0.075365321 0.082455877 0.082455877
## [67] 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877 0.082455877
## [73] 0.082455877 0.082455877 0.089674197 0.096968283 0.104306718 0.104306718
## [79] 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718 0.104306718
## [85] 0.104306718 0.111695071 0.111695071 0.111695071 0.111695071 0.111695071
## [91] 0.111695071 0.126825522 0.126825522 0.126825522 0.126825522 0.126825522
## [97] 0.126825522 0.126825522 0.126825522 0.134984552 0.143227902 0.143227902
## [103] 0.143227902 0.151634521 0.151634521 0.151634521 0.151634521 0.151634521
## [109] 0.151634521 0.151634521 0.160282070 0.160282070 0.168962279 0.168962279
## [115] 0.168962279 0.168962279 0.168962279 0.168962279 0.177720498 0.177720498
## [121] 0.186541816 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [127] 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275 0.195438275
## [133] 0.195438275 0.195438275 0.195438275 0.214054867 0.214054867 0.223603523
## [139] 0.223603523 0.223603523 0.233381313 0.243333679 0.253381794 0.253381794
## [145] 0.253381794 0.253381794 0.253381794 0.263515034 0.263515034 0.263515034
## [151] 0.263515034 0.273710658 0.273710658 0.273710658 0.273710658 0.273710658
## [157] 0.273710658 0.273710658 0.283960762 0.283960762 0.283960762 0.283960762
## [163] 0.283960762 0.283960762 0.283960762 0.283960762 0.294431829 0.305000280
## [169] 0.305000280 0.305000280 0.305000280 0.305000280 0.315826690 0.315826690
## [175] 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499 0.326983499
## [181] 0.326983499 0.338812149 0.338812149 0.338812149 0.363224224 0.363224224
## [187] 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224 0.363224224
## [193] 0.363224224 0.363224224 0.376029567 0.376029567 0.376029567 0.376029567
## [199] 0.376029567 0.376029567 0.376029567 0.389106667 0.389106667 0.402380512
## [205] 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813 0.429762813
## [211] 0.429762813 0.429762813 0.429762813 0.444009105 0.444009105 0.444009105
## [217] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [223] 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284 0.473262284
## [229] 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109 0.489049109
## [235] 0.489049109 0.489049109 0.489049109 0.505600045 0.505600045 0.505600045
## [241] 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045 0.505600045
## [247] 0.505600045 0.543062495 0.543062495 0.543062495 0.543062495 0.543062495
## [253] 0.543062495 0.563339269 0.563339269 0.586028574 0.586028574 0.586028574
## [259] 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574 0.586028574
## [265] 0.586028574 0.612099303 0.612099303 0.612099303 0.612099303 0.639620829
## [271] 0.670502874 0.670502874 0.670502874 0.702589472 0.737090896 0.737090896
## [277] 0.818773352 0.818773352 0.906277419 0.906277419 0.906277419 0.959570305
## [283] 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305 0.959570305
## [289] 0.959570305 0.959570305 1.090085142 1.090085142 1.090085142 1.090085142
## [295] 1.090085142 1.090085142 0.000000000 0.005350627 0.005350627 0.010822944
## [301] 0.010822944 0.010822944 0.010822944 0.016454369 0.016454369 0.016454369
## [307] 0.016454369 0.022164713 0.027889662 0.027889662 0.033763774 0.033763774
## [313] 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774 0.033763774
## [319] 0.033763774 0.033763774 0.033763774 0.033763774 0.039791244 0.039791244
## [325] 0.039791244 0.039791244 0.039791244 0.045954956 0.045954956 0.045954956
## [331] 0.045954956 0.045954956 0.045954956 0.045954956 0.045954956 0.052253850
## [337] 0.052253850 0.052253850 0.052253850 0.058725586 0.058725586 0.058725586
## [343] 0.065226295 0.065226295 0.065226295 0.065226295 0.065226295 0.071800183
## [349] 0.071800183 0.071800183 0.071800183 0.071800183 0.078528422 0.078528422
## [355] 0.085296520 0.085296520 0.085296520 0.085296520 0.092181968 0.099100984
## [361] 0.106076035 0.106076035 0.113136184 0.113136184 0.120399425 0.127705587
## [367] 0.127705587 0.127705587 0.135076428 0.142522708 0.142522708 0.149993842
## [373] 0.157531961 0.157531961 0.157531961 0.157531961 0.165150610 0.165150610
## [379] 0.172846201 0.180620008 0.180620008 0.188445918 0.188445918 0.188445918
## [385] 0.188445918 0.188445918 0.196368910 0.196368910 0.196368910 0.196368910
## [391] 0.204695686 0.204695686 0.204695686 0.213256293 0.213256293 0.213256293
## [397] 0.213256293 0.213256293 0.221921101 0.221921101 0.230668576 0.239491574
## [403] 0.239491574 0.248425299 0.248425299 0.248425299 0.248425299 0.248425299
## [409] 0.257512175 0.266684841 0.266684841 0.266684841 0.266684841 0.266684841
## [415] 0.266684841 0.275887140 0.275887140 0.285301563 0.285301563 0.294762418
## [421] 0.304319540 0.304319540 0.304319540 0.304319540 0.304319540 0.304319540
## [427] 0.304319540 0.304319540 0.304319540 0.314131316 0.314131316 0.324057279
## [433] 0.334041522 0.334041522 0.344145476 0.344145476 0.344145476 0.344145476
## [439] 0.344145476 0.344145476 0.344145476 0.354489364 0.354489364 0.364955963
## [445] 0.364955963 0.375468287 0.375468287 0.375468287 0.375468287 0.375468287
## [451] 0.386067530 0.396810622 0.396810622 0.396810622 0.418872475 0.418872475
## [457] 0.430217353 0.430217353 0.430217353 0.441771494 0.453445406 0.453445406
## [463] 0.465319253 0.465319253 0.477384731 0.489634257 0.502028814 0.514632492
## [469] 0.527340450 0.540164985 0.540164985 0.540164985 0.540164985 0.553304597
## [475] 0.553304597 0.566657491 0.580311979 0.580311979 0.580311979 0.594305960
## [481] 0.609012152 0.623904225 0.639143157 0.639143157 0.639143157 0.654807088
## [487] 0.670718435 0.670718435 0.687301167 0.704071054 0.704071054 0.704071054
## [493] 0.721157987 0.721157987 0.738561528 0.756346779 0.756346779 0.756346779
## [499] 0.774488197 0.774488197 0.774488197 0.774488197 0.774488197 0.793145557
## [505] 0.793145557 0.812226345 0.812226345 0.812226345 0.831644405 0.851361126
## [511] 0.871448014 0.871448014 0.871448014 0.891810367 0.891810367 0.913268483
## [517] 0.913268483 0.913268483 0.913268483 0.913268483 0.936798339 0.936798339
## [523] 0.936798339 0.936798339 0.936798339 0.961841394 0.987712859 0.987712859
## [529] 1.013960311 1.040953046 1.040953046 1.040953046 1.069941379 1.099543534
## [535] 1.099543534 1.099543534 1.099543534 1.099543534 1.129531930 1.129531930
## [541] 1.160547040 1.192692166 1.228912605 1.228912605 1.228912605 1.228912605
## [547] 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605 1.228912605
## [553] 1.228912605 1.270816176 1.314046197 1.314046197 1.359029219 1.404934747
## [559] 1.404934747 1.404934747 1.404934747 1.404934747 1.458863406 1.458863406
## [565] 1.458863406 1.458863406 1.520684154 1.520684154 1.590254427 1.590254427
## [571] 1.667260347 1.667260347 1.747291122 1.747291122 1.747291122 1.747291122
## [577] 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122 1.747291122
## [583] 1.747291122 1.747291122 1.747291122 1.926625374 1.926625374 1.926625374
## [589] 1.926625374 1.926625374 2.281783691 2.281783691 0.009633169 0.009633169
## [595] 0.014551755 0.014551755 0.014551755 0.019546218 0.019546218 0.024599344
## [601] 0.029678509 0.034864290 0.040107303 0.045413686 0.045413686 0.050750288
## [607] 0.050750288 0.061572580 0.067117647 0.072742518 0.078415489 0.078415489
## [613] 0.084121378 0.089991681 0.095877420 0.101853535 0.107873484 0.107873484
## [619] 0.107873484 0.113924256 0.119989528 0.119989528 0.126231259 0.126231259
## [625] 0.132564626 0.132564626 0.138970374 0.138970374 0.138970374 0.138970374
## [631] 0.158792804 0.158792804 0.165618521 0.165618521 0.172484282 0.172484282
## [637] 0.186450350 0.193667471 0.200947081 0.208250442 0.223100168 0.230731258
## [643] 0.230731258 0.230731258 0.238451311 0.246241191 0.246241191 0.246241191
## [649] 0.246241191 0.254091780 0.254091780 0.261998834 0.269933185 0.269933185
## [655] 0.295740183 0.295740183 0.304607658 0.304607658 0.304607658 0.313541365
## [661] 0.313541365 0.322590148 0.331761267 0.341155661 0.341155661 0.350635123
## [667] 0.350635123 0.360361864 0.360361864 0.370472540 0.380732678 0.391154668
## [673] 0.391154668 0.401655837 0.401655837 0.401655837 0.412274302 0.412274302
## [679] 0.422970130 0.422970130 0.422970130 0.433912349 0.433912349 0.433912349
## [685] 0.433912349 0.445124076 0.445124076 0.445124076 0.456660219 0.456660219
## [691] 0.456660219 0.468424656 0.480317909 0.480317909 0.480317909 0.492387762
## [697] 0.492387762 0.492387762 0.504593246 0.504593246 0.516909136 0.516909136
## [703] 0.529347066 0.529347066 0.529347066 0.529347066 0.541876644 0.554515396
## [709] 0.567272570 0.580175913 0.593186674 0.593186674 0.606344153 0.606344153
## [715] 0.619675763 0.619675763 0.619675763 0.633344103 0.647161102 0.647161102
## [721] 0.647161102 0.647161102 0.661946165 0.677209373 0.677209373 0.677209373
## [727] 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700 0.692791700
## [733] 0.708979609 0.708979609 0.725379668 0.725379668 0.742673869 0.742673869
## [739] 0.760665334 0.760665334 0.760665334 0.760665334 0.779916835 0.799534319
## [745] 0.819325249 0.839297067 0.839297067 0.839297067 0.859999536 0.859999536
## [751] 0.859999536 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673
## [757] 0.881788673 0.881788673 0.881788673 0.881788673 0.881788673 0.905686340
## [763] 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340 0.905686340
## [769] 0.930593182 0.930593182 0.956024835 0.956024835 0.956024835 0.956024835
## [775] 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461 0.982860461
## [781] 0.982860461 0.982860461 0.982860461 1.011138349 1.011138349 1.011138349
## [787] 1.040995345 1.072147856 1.072147856 1.104113664 1.104113664 1.104113664
## [793] 1.136775819 1.171345121 1.171345121 1.207040278 1.207040278 1.207040278
## [799] 1.207040278 1.207040278 1.246907275 1.246907275 1.287930715 1.330411123
## [805] 1.330411123 1.330411123 1.330411123 1.374604413 1.420519943 1.420519943
## [811] 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943 1.420519943
## [817] 1.420519943 1.420519943 1.470281401 1.521709021 1.521709021 1.521709021
## [823] 1.521709021 1.574890670 1.574890670 1.574890670 1.574890670 1.574890670
## [829] 1.574890670 1.574890670 1.629630555 1.685791093 1.685791093 1.685791093
## [835] 1.685791093 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [841] 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642 1.750311642
## [847] 1.832400858 1.832400858 1.917615102 1.917615102 1.917615102 1.917615102
## [853] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [859] 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102 1.917615102
## [865] 1.917615102 2.029441753 2.029441753 2.029441753 2.029441753 2.029441753
## [871] 2.029441753 2.029441753 2.029441753 2.029441753 2.194031226 2.194031226
## [877] 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226 2.194031226
## [883] 2.530841239 2.925483659 2.925483659 2.925483659 2.925483659 2.925483659
res4<-residuals(Cox.fit4,type="deviance")
hist(res4)
<<<<<<< HEAD
boxplot(res4~Cdata$Death_Censoring)
boxplot(res4~Cdata$Death_Censoring)
# Looking at Protrombin
min(Cdata[,8])# Lowest 12
## [1] 12
max(Cdata[,8])# Highest 135
## [1] 135
Prothrombin.groups=cut(Cdata$Prothrombin,breaks=c(-Inf,quantile(Cdata$Prothrombin ,prob=c(0.33,0.66)),Inf)) #devides the Age into groups
Cox.fit5 <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ strata(Prothrombin.groups), data = Cdata)
summary(Cox.fit5)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + strata(Prothrombin.groups), data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.100662 1.105902 0.118140 0.852 0.39418
## factor(Sex)1 0.475374 1.608615 0.126927 3.745 0.00018 ***
## factor(Ascites)1 0.460651 1.585105 0.178184 2.585 0.00973 **
## factor(Ascites)2 0.952334 2.591752 0.184248 5.169 2.36e-07 ***
## Age 0.048430 1.049622 0.006825 7.096 1.28e-12 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.106 0.9042 0.8773 1.394
## factor(Sex)1 1.609 0.6217 1.2543 2.063
## factor(Ascites)1 1.585 0.6309 1.1179 2.248
## factor(Ascites)2 2.592 0.3858 1.8062 3.719
## Age 1.050 0.9527 1.0357 1.064
##
## Concordance= 0.672 (se = 0.018 )
## Likelihood ratio test= 91.68 on 5 df, p=<2e-16
## Wald test = 90.79 on 5 df, p=<2e-16
## Score (logrank) test = 93.72 on 5 df, p=<2e-16
source("andersen.R")
test5<-survfit(Cox.fit5)
test5
## Call: survfit(formula = Cox.fit5)
##
## n events median 0.95LCL 0.95UCL
## (-Inf,58] 171 127 1169 883 1582
## (58,80] 165 96 1698 1203 2195
## (80, Inf] 152 69 2364 1976 3132
andersen.plot(test5) # must have strata, plots the cumilated hazard function for the levels of the strata variable
<<<<<<< HEAD
## [1] 0.003815209 0.007683007 0.011584233 0.015524640 0.019505384 0.023538256
## [7] 0.023538256 0.027640426 0.035938955 0.035938955 0.040147080 0.044419856
## [13] 0.044419856 0.044419856 0.044419856 0.044419856 0.048804220 0.062254913
## [19] 0.066857759 0.076224418 0.076224418 0.076224418 0.076224418 0.080982967
## [25] 0.085819661 0.090716044 0.095649157 0.095649157 0.100605383 0.100605383
## [31] 0.105595367 0.105595367 0.110675694 0.115868619 0.115868619 0.115868619
## [37] 0.121148141 0.121148141 0.121148141 0.121148141 0.121148141 0.131893456
## [43] 0.131893456 0.137421743 0.142988787 0.142988787 0.148646614 0.148646614
## [49] 0.160228415 0.160228415 0.166295953 0.178610032 0.178610032 0.191138887
## [55] 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449
## [61] 0.204138206 0.210707735 0.217318317 0.224215818 0.238147577 0.238147577
## [67] 0.245219522 0.252306963 0.259472322 0.259472322 0.274128676 0.274128676
## [73] 0.274128676 0.281546394 0.281546394 0.289018618 0.289018618 0.296531974
## [79] 0.304353874 0.312209872 0.312209872 0.320138856 0.320138856 0.320138856
## [85] 0.328159980 0.336243939 0.336243939 0.344381675 0.344381675 0.352598926
## [91] 0.352598926 0.360910937 0.369392228 0.378042605 0.378042605 0.378042605
## [97] 0.378042605 0.386845320 0.386845320 0.395828950 0.395828950 0.395828950
## [103] 0.395828950 0.395828950 0.404948309 0.414138672 0.414138672 0.414138672
## [109] 0.414138672 0.414138672 0.423426979 0.423426979 0.423426979 0.432806013
## [115] 0.432806013 0.442290826 0.442290826 0.451875011 0.451875011 0.461510440
## [121] 0.461510440 0.471258993 0.471258993 0.471258993 0.481150204 0.481150204
## [127] 0.481150204 0.481150204 0.491080246 0.491080246 0.491080246 0.491080246
## [133] 0.491080246 0.501355641 0.511878399 0.511878399 0.511878399 0.522589991
## [139] 0.522589991 0.533493357 0.533493357 0.544539943 0.544539943 0.544539943
## [145] 0.544539943 0.544539943 0.544539943 0.544539943 0.555981864 0.567486954
## [151] 0.567486954 0.579419841 0.579419841 0.579419841 0.579419841 0.579419841
## [157] 0.579419841 0.592453938 0.605747671 0.605747671 0.619209067 0.633357923
## [163] 0.633357923 0.633357923 0.633357923 0.633357923 0.633357923 0.648497735
## [169] 0.663809559 0.663809559 0.663809559 0.663809559 0.663809559 0.679715552
## [175] 0.695849750 0.712319855 0.729088039 0.729088039 0.729088039 0.729088039
## [181] 0.729088039 0.746646274 0.764560852 0.764560852 0.764560852 0.783263260
## [187] 0.802337188 0.802337188 0.822432633 0.842613180 0.842613180 0.863519292
## [193] 0.863519292 0.863519292 0.884708578 0.884708578 0.884708578 0.884708578
## [199] 0.906698043 0.906698043 0.928838016 0.928838016 0.928838016 0.928838016
## [205] 0.928838016 0.951641911 0.951641911 0.951641911 0.951641911 0.951641911
## [211] 0.951641911 0.951641911 0.951641911 0.975097788 0.999170914 0.999170914
## [217] 1.023722913 1.048726373 1.074879725 1.074879725 1.102723976 1.102723976
## [223] 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976
## [229] 1.132454776 1.132454776 1.132454776 1.132454776 1.163457833 1.196856365
## [235] 1.196856365 1.196856365 1.231339476 1.231339476 1.271278329 1.271278329
## [241] 1.271278329 1.315576550 1.315576550 1.360742931 1.360742931 1.360742931
## [247] 1.360742931 1.408744019 1.408744019 1.408744019 1.408744019 1.408744019
## [253] 1.408744019 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [259] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [265] 1.461889660 1.461889660 1.524568595 1.524568595 1.594403830 1.594403830
## [271] 1.668862962 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [277] 1.751026186 1.751026186 1.842255469 1.842255469 1.842255469 1.842255469
## [283] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469
## [289] 1.842255469 1.842255469 1.842255469 1.842255469 2.124490838 2.124490838
## [295] 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161
## [301] 2.496124161 2.496124161 2.496124161 0.000000000 0.000000000 0.000000000
## [307] 0.000000000 0.000000000 0.000000000 0.000000000 0.005214536 0.005214536
## [313] 0.005214536 0.005214536 0.005214536 0.010551942 0.010551942 0.015947912
## [319] 0.015947912 0.015947912 0.015947912 0.015947912 0.015947912 0.021407242
## [325] 0.021407242 0.021407242 0.021407242 0.021407242 0.021407242 0.021407242
## [331] 0.021407242 0.021407242 0.026971836 0.026971836 0.026971836 0.026971836
## [337] 0.026971836 0.032611990 0.032611990 0.032611990 0.038280941 0.038280941
## [343] 0.038280941 0.038280941 0.044119420 0.050049846 0.050049846 0.050049846
## [349] 0.050049846 0.050049846 0.050049846 0.050049846 0.056090678 0.056090678
## [355] 0.056090678 0.062210413 0.062210413 0.062210413 0.062210413 0.068421665
## [361] 0.074704142 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806
## [367] 0.081095806 0.094314691 0.101097016 0.101097016 0.101097016 0.101097016
## [373] 0.101097016 0.101097016 0.107965714 0.115044736 0.115044736 0.122190020
## [379] 0.122190020 0.122190020 0.129556960 0.129556960 0.129556960 0.136977660
## [385] 0.136977660 0.144630677 0.152345180 0.152345180 0.152345180 0.160176278
## [391] 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278
## [397] 0.160176278 0.168197928 0.168197928 0.168197928 0.168197928 0.168197928
## [403] 0.168197928 0.168197928 0.168197928 0.176614537 0.176614537 0.176614537
## [409] 0.176614537 0.185176955 0.193829682 0.193829682 0.202671703 0.211594856
## [415] 0.220626646 0.220626646 0.220626646 0.229793501 0.229793501 0.239077829
## [421] 0.248444048 0.257932939 0.257932939 0.267521287 0.267521287 0.277197004
## [427] 0.286917451 0.286917451 0.286917451 0.296920891 0.307216111 0.307216111
## [433] 0.317632823 0.317632823 0.317632823 0.317632823 0.317632823 0.317632823
## [439] 0.328847214 0.340199006 0.340199006 0.351912350 0.351912350 0.351912350
## [445] 0.351912350 0.351912350 0.363932510 0.363932510 0.376103680 0.388411924
## [451] 0.388411924 0.388411924 0.388411924 0.388411924 0.388411924 0.400827235
## [457] 0.413388661 0.426093822 0.426093822 0.438880811 0.438880811 0.438880811
## [463] 0.451909153 0.451909153 0.465106583 0.478494715 0.492428623 0.492428623
## [469] 0.492428623 0.492428623 0.492428623 0.492428623 0.492428623 0.506710060
## [475] 0.521152134 0.521152134 0.521152134 0.535918656 0.535918656 0.535918656
## [481] 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093
## [487] 0.566225029 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [493] 0.581917832 0.581917832 0.581917832 0.581917832 0.598575995 0.598575995
## [499] 0.615540328 0.633502999 0.651685645 0.651685645 0.670216074 0.670216074
## [505] 0.689056961 0.708270080 0.728416190 0.748931218 0.748931218 0.769827448
## [511] 0.791444671 0.813230253 0.813230253 0.835981763 0.859337551 0.859337551
## [517] 0.859337551 0.859337551 0.883790398 0.883790398 0.883790398 0.883790398
## [523] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [529] 0.910092808 0.937699114 0.966355040 0.966355040 0.995929592 1.026020500
## [535] 1.057173383 1.057173383 1.057173383 1.057173383 1.091579070 1.091579070
## [541] 1.127017733 1.127017733 1.127017733 1.163673474 1.163673474 1.163673474
## [547] 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474
## [553] 1.163673474 1.205390594 1.205390594 1.205390594 1.205390594 1.248789143
## [559] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [565] 1.248789143 1.302049650 1.302049650 1.302049650 1.302049650 1.361380792
## [571] 1.361380792 1.361380792 1.423171756 1.423171756 1.423171756 1.423171756
## [577] 1.423171756 1.492959677 1.565134027 1.641566599 1.641566599 1.641566599
## [583] 1.641566599 1.641566599 1.641566599 1.760594965 1.760594965 1.760594965
## [589] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [595] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.995646102
## [601] 1.995646102 1.995646102 1.995646102 1.995646102 1.995646102 1.995646102
## [607] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [613] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [619] 0.007077094 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278
## [625] 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278
## [631] 0.021477375 0.021477375 0.021477375 0.021477375 0.021477375 0.021477375
## [637] 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193
## [643] 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193 0.043683254
## [649] 0.043683254 0.043683254 0.043683254 0.043683254 0.043683254 0.051238198
## [655] 0.051238198 0.058848646 0.058848646 0.066482062 0.066482062 0.066482062
## [661] 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062
## [667] 0.066482062 0.066482062 0.082161595 0.082161595 0.082161595 0.082161595
## [673] 0.082161595 0.090312041 0.090312041 0.107058768 0.107058768 0.107058768
## [679] 0.107058768 0.107058768 0.115642909 0.115642909 0.115642909 0.124279460
## [685] 0.124279460 0.124279460 0.124279460 0.124279460 0.124279460 0.132971217
## [691] 0.132971217 0.132971217 0.132971217 0.141833373 0.141833373 0.141833373
## [697] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [703] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [709] 0.141833373 0.141833373 0.141833373 0.151197061 0.151197061 0.151197061
## [715] 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061
## [721] 0.160686642 0.160686642 0.170396775 0.170396775 0.180213133 0.190108901
## [727] 0.190108901 0.210270329 0.210270329 0.210270329 0.210270329 0.210270329
## [733] 0.220576216 0.220576216 0.230965174 0.230965174 0.230965174 0.230965174
## [739] 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174
## [745] 0.241881805 0.241881805 0.241881805 0.241881805 0.241881805 0.241881805
## [751] 0.241881805 0.241881805 0.252853192 0.252853192 0.252853192 0.252853192
## [757] 0.252853192 0.252853192 0.252853192 0.252853192 0.252853192 0.264161947
## [763] 0.275576224 0.310613892 0.310613892 0.310613892 0.310613892 0.322632099
## [769] 0.322632099 0.322632099 0.322632099 0.322632099 0.322632099 0.322632099
## [775] 0.322632099 0.322632099 0.334704748 0.334704748 0.334704748 0.334704748
## [781] 0.334704748 0.334704748 0.347024049 0.372137832 0.372137832 0.372137832
## [787] 0.372137832 0.385220199 0.385220199 0.385220199 0.398660095 0.398660095
## [793] 0.398660095 0.398660095 0.412484044 0.426505935 0.426505935 0.426505935
## [799] 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935
## [805] 0.426505935 0.456128884 0.471685153 0.471685153 0.471685153 0.471685153
## [811] 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692
## [817] 0.487657692 0.504636398 0.504636398 0.504636398 0.540398037 0.540398037
## [823] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442
## [829] 0.597092442 0.597092442 0.597092442 0.597092442 0.617805328 0.617805328
## [835] 0.617805328 0.617805328 0.640376113 0.640376113 0.664790137 0.715491227
## [841] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [847] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [853] 0.715491227 0.715491227 0.715491227 0.715491227 0.746629486 0.746629486
## [859] 0.746629486 0.746629486 0.746629486 0.746629486 0.780446351 0.780446351
## [865] 0.851972559 0.851972559 0.934636454 0.934636454 0.934636454 0.934636454
## [871] 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454
## [877] 0.934636454 0.934636454 0.934636454 0.986212450 0.986212450 1.103396433
## [883] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [889] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [895] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.238566987
## [901] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
## [907] 1.390361656 1.390361656 1.390361656
## [1] 0.003815209 0.007683007 0.011584233 0.015524640 0.019505384 0.023538256
## [7] 0.023538256 0.023538256 0.027640426 0.035938955 0.035938955 0.040147080
## [13] 0.044419856 0.044419856 0.044419856 0.044419856 0.044419856 0.048804220
## [19] 0.062254913 0.066857759 0.076224418 0.076224418 0.080982967 0.080982967
## [25] 0.085819661 0.090716044 0.095649157 0.100605383 0.100605383 0.100605383
## [31] 0.105595367 0.110675694 0.115868619 0.115868619 0.121148141 0.121148141
## [37] 0.131893456 0.137421743 0.142988787 0.142988787 0.148646614 0.148646614
## [43] 0.160228415 0.160228415 0.166295953 0.178610032 0.191138887 0.197592449
## [49] 0.197592449 0.197592449 0.197592449 0.204138206 0.210707735 0.210707735
## [55] 0.210707735 0.210707735 0.217318317 0.224215818 0.238147577 0.245219522
## [61] 0.252306963 0.252306963 0.259472322 0.259472322 0.259472322 0.274128676
## [67] 0.274128676 0.281546394 0.281546394 0.289018618 0.289018618 0.296531974
## [73] 0.304353874 0.312209872 0.320138856 0.320138856 0.320138856 0.328159980
## [79] 0.336243939 0.344381675 0.352598926 0.352598926 0.360910937 0.369392228
## [85] 0.378042605 0.386845320 0.395828950 0.395828950 0.404948309 0.404948309
## [91] 0.404948309 0.414138672 0.414138672 0.414138672 0.423426979 0.432806013
## [97] 0.432806013 0.442290826 0.442290826 0.442290826 0.442290826 0.442290826
## [103] 0.451875011 0.451875011 0.451875011 0.461510440 0.461510440 0.461510440
## [109] 0.461510440 0.471258993 0.481150204 0.481150204 0.481150204 0.491080246
## [115] 0.491080246 0.491080246 0.501355641 0.511878399 0.511878399 0.522589991
## [121] 0.522589991 0.533493357 0.544539943 0.544539943 0.544539943 0.544539943
## [127] 0.555981864 0.567486954 0.567486954 0.579419841 0.579419841 0.579419841
## [133] 0.579419841 0.579419841 0.579419841 0.579419841 0.592453938 0.605747671
## [139] 0.619209067 0.619209067 0.619209067 0.633357923 0.633357923 0.633357923
## [145] 0.648497735 0.663809559 0.663809559 0.663809559 0.679715552 0.679715552
## [151] 0.695849750 0.712319855 0.712319855 0.712319855 0.729088039 0.729088039
## [157] 0.729088039 0.729088039 0.729088039 0.729088039 0.729088039 0.729088039
## [163] 0.746646274 0.764560852 0.764560852 0.764560852 0.783263260 0.802337188
## [169] 0.802337188 0.822432633 0.822432633 0.822432633 0.842613180 0.842613180
## [175] 0.863519292 0.863519292 0.884708578 0.906698043 0.906698043 0.906698043
## [181] 0.906698043 0.906698043 0.928838016 0.928838016 0.951641911 0.951641911
## [187] 0.951641911 0.951641911 0.975097788 0.975097788 0.975097788 0.975097788
## [193] 0.999170914 0.999170914 0.999170914 0.999170914 1.023722913 1.048726373
## [199] 1.074879725 1.074879725 1.102723976 1.102723976 1.102723976 1.102723976
## [205] 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976
## [211] 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776
## [217] 1.132454776 1.163457833 1.163457833 1.163457833 1.196856365 1.196856365
## [223] 1.196856365 1.231339476 1.271278329 1.271278329 1.315576550 1.360742931
## [229] 1.360742931 1.360742931 1.360742931 1.360742931 1.360742931 1.360742931
## [235] 1.408744019 1.408744019 1.408744019 1.408744019 1.408744019 1.408744019
## [241] 1.408744019 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [247] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [253] 1.461889660 1.524568595 1.524568595 1.594403830 1.594403830 1.668862962
## [259] 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [265] 1.751026186 1.751026186 1.751026186 1.751026186 1.842255469 1.842255469
## [271] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469
## [277] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 2.124490838
## [283] 2.124490838 2.124490838 2.124490838 2.496124161 2.496124161 2.496124161
## [289] 2.496124161 2.496124161 2.496124161 2.496124161 0.000000000 0.000000000
## [295] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [301] 0.005214536 0.005214536 0.005214536 0.005214536 0.005214536 0.010551942
## [307] 0.010551942 0.010551942 0.015947912 0.015947912 0.015947912 0.015947912
## [313] 0.015947912 0.021407242 0.021407242 0.021407242 0.021407242 0.021407242
## [319] 0.021407242 0.021407242 0.026971836 0.026971836 0.026971836 0.026971836
## [325] 0.026971836 0.032611990 0.032611990 0.038280941 0.044119420 0.050049846
## [331] 0.050049846 0.050049846 0.050049846 0.050049846 0.050049846 0.050049846
## [337] 0.056090678 0.056090678 0.062210413 0.062210413 0.081095806 0.081095806
## [343] 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806
## [349] 0.081095806 0.081095806 0.094314691 0.101097016 0.101097016 0.101097016
## [355] 0.101097016 0.101097016 0.101097016 0.101097016 0.115044736 0.115044736
## [361] 0.122190020 0.122190020 0.122190020 0.129556960 0.129556960 0.129556960
## [367] 0.136977660 0.144630677 0.152345180 0.152345180 0.152345180 0.160176278
## [373] 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278 0.168197928
## [379] 0.168197928 0.168197928 0.176614537 0.176614537 0.176614537 0.176614537
## [385] 0.176614537 0.202671703 0.211594856 0.220626646 0.220626646 0.229793501
## [391] 0.229793501 0.229793501 0.229793501 0.239077829 0.248444048 0.248444048
## [397] 0.257932939 0.257932939 0.267521287 0.267521287 0.267521287 0.267521287
## [403] 0.286917451 0.286917451 0.307216111 0.307216111 0.317632823 0.317632823
## [409] 0.317632823 0.317632823 0.317632823 0.340199006 0.340199006 0.351912350
## [415] 0.351912350 0.363932510 0.376103680 0.388411924 0.388411924 0.388411924
## [421] 0.388411924 0.388411924 0.426093822 0.426093822 0.426093822 0.438880811
## [427] 0.438880811 0.438880811 0.438880811 0.438880811 0.451909153 0.451909153
## [433] 0.451909153 0.465106583 0.492428623 0.492428623 0.492428623 0.492428623
## [439] 0.492428623 0.492428623 0.521152134 0.521152134 0.535918656 0.535918656
## [445] 0.535918656 0.535918656 0.535918656 0.535918656 0.535918656 0.550949093
## [451] 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093 0.550949093
## [457] 0.566225029 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [463] 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [469] 0.598575995 0.651685645 0.651685645 0.651685645 0.670216074 0.670216074
## [475] 0.670216074 0.728416190 0.748931218 0.813230253 0.813230253 0.835981763
## [481] 0.859337551 0.859337551 0.859337551 0.859337551 0.859337551 0.883790398
## [487] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [493] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [499] 0.910092808 0.910092808 0.966355040 0.966355040 0.966355040 0.995929592
## [505] 0.995929592 0.995929592 1.057173383 1.057173383 1.057173383 1.057173383
## [511] 1.057173383 1.057173383 1.057173383 1.057173383 1.091579070 1.091579070
## [517] 1.127017733 1.127017733 1.163673474 1.163673474 1.163673474 1.163673474
## [523] 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474
## [529] 1.163673474 1.163673474 1.163673474 1.205390594 1.205390594 1.205390594
## [535] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [541] 1.248789143 1.248789143 1.248789143 1.248789143 1.302049650 1.361380792
## [547] 1.361380792 1.361380792 1.423171756 1.423171756 1.423171756 1.423171756
## [553] 1.423171756 1.423171756 1.492959677 1.492959677 1.492959677 1.641566599
## [559] 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599
## [565] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [571] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [577] 1.760594965 1.760594965 1.760594965 1.760594965 1.995646102 1.995646102
## [583] 1.995646102 1.995646102 0.000000000 0.000000000 0.000000000 0.000000000
## [589] 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000 0.000000000
## [595] 0.000000000 0.000000000 0.000000000 0.007077094 0.007077094 0.014216278
## [601] 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278 0.014216278
## [607] 0.014216278 0.021477375 0.021477375 0.021477375 0.021477375 0.021477375
## [613] 0.028820277 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193
## [619] 0.036190193 0.036190193 0.043683254 0.043683254 0.043683254 0.043683254
## [625] 0.043683254 0.051238198 0.051238198 0.058848646 0.058848646 0.066482062
## [631] 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062
## [637] 0.066482062 0.066482062 0.074283328 0.082161595 0.082161595 0.082161595
## [643] 0.082161595 0.082161595 0.090312041 0.090312041 0.090312041 0.098605212
## [649] 0.107058768 0.107058768 0.107058768 0.107058768 0.115642909 0.115642909
## [655] 0.124279460 0.124279460 0.124279460 0.124279460 0.124279460 0.132971217
## [661] 0.132971217 0.132971217 0.132971217 0.141833373 0.141833373 0.141833373
## [667] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [673] 0.141833373 0.141833373 0.151197061 0.151197061 0.151197061 0.151197061
## [679] 0.151197061 0.151197061 0.160686642 0.160686642 0.160686642 0.160686642
## [685] 0.170396775 0.170396775 0.170396775 0.180213133 0.190108901 0.190108901
## [691] 0.200148605 0.210270329 0.210270329 0.210270329 0.210270329 0.220576216
## [697] 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174
## [703] 0.230965174 0.230965174 0.241881805 0.241881805 0.241881805 0.241881805
## [709] 0.252853192 0.252853192 0.252853192 0.252853192 0.252853192 0.252853192
## [715] 0.264161947 0.264161947 0.275576224 0.287173209 0.298818473 0.310613892
## [721] 0.310613892 0.310613892 0.310613892 0.310613892 0.322632099 0.322632099
## [727] 0.322632099 0.322632099 0.322632099 0.322632099 0.322632099 0.334704748
## [733] 0.334704748 0.334704748 0.334704748 0.334704748 0.334704748 0.347024049
## [739] 0.347024049 0.359487200 0.372137832 0.372137832 0.372137832 0.372137832
## [745] 0.372137832 0.385220199 0.385220199 0.385220199 0.398660095 0.398660095
## [751] 0.398660095 0.398660095 0.412484044 0.412484044 0.412484044 0.426505935
## [757] 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935
## [763] 0.441039991 0.456128884 0.456128884 0.471685153 0.471685153 0.487657692
## [769] 0.487657692 0.487657692 0.487657692 0.504636398 0.504636398 0.504636398
## [775] 0.522274168 0.540398037 0.540398037 0.558903304 0.577733813 0.597092442
## [781] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442
## [787] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.617805328
## [793] 0.617805328 0.617805328 0.617805328 0.617805328 0.617805328 0.640376113
## [799] 0.640376113 0.640376113 0.664790137 0.664790137 0.689971460 0.715491227
## [805] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [811] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [817] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [823] 0.746629486 0.746629486 0.746629486 0.746629486 0.780446351 0.815837029
## [829] 0.851972559 0.851972559 0.851972559 0.851972559 0.892199662 0.934636454
## [835] 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454
## [841] 0.934636454 0.934636454 0.934636454 0.934636454 0.986212450 0.986212450
## [847] 1.043419749 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [853] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [859] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [865] 1.103396433 1.103396433 1.238566987 1.238566987 1.390361656 1.390361656
## [871] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
## [1] 0.023538256 0.027640426 0.035938955 0.035938955 0.044419856 0.044419856
## [7] 0.044419856 0.044419856 0.044419856 0.076224418 0.076224418 0.080982967
## [13] 0.095649157 0.100605383 0.100605383 0.100605383 0.105595367 0.110675694
## [19] 0.115868619 0.121148141 0.121148141 0.121148141 0.121148141 0.131893456
## [25] 0.131893456 0.148646614 0.160228415 0.160228415 0.178610032 0.178610032
## [31] 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449 0.197592449
## [37] 0.210707735 0.210707735 0.210707735 0.217318317 0.224215818 0.238147577
## [43] 0.238147577 0.252306963 0.252306963 0.259472322 0.259472322 0.259472322
## [49] 0.274128676 0.274128676 0.274128676 0.281546394 0.289018618 0.289018618
## [55] 0.296531974 0.312209872 0.320138856 0.320138856 0.320138856 0.336243939
## [61] 0.344381675 0.344381675 0.352598926 0.352598926 0.378042605 0.378042605
## [67] 0.378042605 0.386845320 0.395828950 0.395828950 0.395828950 0.395828950
## [73] 0.395828950 0.404948309 0.404948309 0.414138672 0.414138672 0.414138672
## [79] 0.414138672 0.414138672 0.414138672 0.423426979 0.423426979 0.423426979
## [85] 0.432806013 0.432806013 0.442290826 0.442290826 0.442290826 0.442290826
## [91] 0.442290826 0.451875011 0.451875011 0.451875011 0.451875011 0.461510440
## [97] 0.461510440 0.461510440 0.461510440 0.471258993 0.471258993 0.481150204
## [103] 0.481150204 0.481150204 0.481150204 0.481150204 0.491080246 0.491080246
## [109] 0.491080246 0.491080246 0.501355641 0.511878399 0.511878399 0.511878399
## [115] 0.522589991 0.522589991 0.533493357 0.544539943 0.544539943 0.544539943
## [121] 0.544539943 0.544539943 0.579419841 0.579419841 0.579419841 0.579419841
## [127] 0.579419841 0.579419841 0.579419841 0.579419841 0.579419841 0.605747671
## [133] 0.619209067 0.619209067 0.633357923 0.633357923 0.633357923 0.633357923
## [139] 0.663809559 0.663809559 0.663809559 0.663809559 0.679715552 0.695849750
## [145] 0.712319855 0.712319855 0.729088039 0.729088039 0.729088039 0.729088039
## [151] 0.729088039 0.729088039 0.764560852 0.764560852 0.764560852 0.764560852
## [157] 0.783263260 0.802337188 0.802337188 0.822432633 0.822432633 0.863519292
## [163] 0.863519292 0.863519292 0.884708578 0.884708578 0.884708578 0.906698043
## [169] 0.906698043 0.906698043 0.906698043 0.906698043 0.928838016 0.928838016
## [175] 0.928838016 0.928838016 0.928838016 0.951641911 0.951641911 0.951641911
## [181] 0.951641911 0.951641911 0.951641911 0.951641911 0.951641911 0.951641911
## [187] 0.951641911 0.975097788 0.975097788 0.975097788 0.999170914 0.999170914
## [193] 0.999170914 0.999170914 1.102723976 1.102723976 1.102723976 1.102723976
## [199] 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976 1.102723976
## [205] 1.102723976 1.102723976 1.132454776 1.132454776 1.132454776 1.132454776
## [211] 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776 1.132454776
## [217] 1.163457833 1.163457833 1.196856365 1.196856365 1.196856365 1.196856365
## [223] 1.231339476 1.271278329 1.315576550 1.360742931 1.360742931 1.360742931
## [229] 1.360742931 1.360742931 1.360742931 1.360742931 1.408744019 1.408744019
## [235] 1.408744019 1.408744019 1.408744019 1.461889660 1.461889660 1.461889660
## [241] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [247] 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660 1.461889660
## [253] 1.461889660 1.461889660 1.461889660 1.461889660 1.524568595 1.594403830
## [259] 1.594403830 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [265] 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186 1.751026186
## [271] 1.751026186 1.751026186 1.842255469 1.842255469 1.842255469 1.842255469
## [277] 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469 1.842255469
## [283] 1.842255469 2.124490838 2.124490838 2.496124161 2.496124161 2.496124161
## [289] 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161 2.496124161
## [295] 2.496124161 0.000000000 0.005214536 0.005214536 0.005214536 0.010551942
## [301] 0.010551942 0.010551942 0.015947912 0.015947912 0.021407242 0.021407242
## [307] 0.021407242 0.021407242 0.026971836 0.026971836 0.026971836 0.026971836
## [313] 0.026971836 0.032611990 0.038280941 0.038280941 0.038280941 0.038280941
## [319] 0.044119420 0.050049846 0.050049846 0.050049846 0.056090678 0.056090678
## [325] 0.062210413 0.062210413 0.068421665 0.074704142 0.081095806 0.081095806
## [331] 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806 0.081095806
## [337] 0.094314691 0.101097016 0.101097016 0.101097016 0.101097016 0.101097016
## [343] 0.101097016 0.107965714 0.115044736 0.115044736 0.122190020 0.122190020
## [349] 0.122190020 0.129556960 0.136977660 0.144630677 0.144630677 0.152345180
## [355] 0.160176278 0.160176278 0.160176278 0.160176278 0.160176278 0.168197928
## [361] 0.168197928 0.168197928 0.168197928 0.168197928 0.168197928 0.168197928
## [367] 0.176614537 0.176614537 0.176614537 0.176614537 0.176614537 0.185176955
## [373] 0.193829682 0.193829682 0.202671703 0.202671703 0.211594856 0.220626646
## [379] 0.220626646 0.220626646 0.229793501 0.229793501 0.229793501 0.229793501
## [385] 0.239077829 0.239077829 0.248444048 0.248444048 0.257932939 0.257932939
## [391] 0.267521287 0.267521287 0.267521287 0.267521287 0.277197004 0.286917451
## [397] 0.286917451 0.286917451 0.296920891 0.307216111 0.307216111 0.317632823
## [403] 0.317632823 0.317632823 0.317632823 0.317632823 0.317632823 0.328847214
## [409] 0.340199006 0.340199006 0.351912350 0.351912350 0.351912350 0.363932510
## [415] 0.376103680 0.376103680 0.388411924 0.400827235 0.413388661 0.426093822
## [421] 0.426093822 0.426093822 0.438880811 0.438880811 0.438880811 0.438880811
## [427] 0.451909153 0.451909153 0.451909153 0.465106583 0.478494715 0.492428623
## [433] 0.492428623 0.492428623 0.506710060 0.521152134 0.521152134 0.521152134
## [439] 0.535918656 0.535918656 0.535918656 0.535918656 0.535918656 0.550949093
## [445] 0.550949093 0.550949093 0.550949093 0.566225029 0.566225029 0.581917832
## [451] 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832 0.581917832
## [457] 0.581917832 0.581917832 0.598575995 0.615540328 0.633502999 0.651685645
## [463] 0.651685645 0.651685645 0.670216074 0.670216074 0.670216074 0.689056961
## [469] 0.708270080 0.728416190 0.728416190 0.748931218 0.769827448 0.791444671
## [475] 0.813230253 0.813230253 0.813230253 0.813230253 0.835981763 0.835981763
## [481] 0.859337551 0.859337551 0.859337551 0.859337551 0.859337551 0.883790398
## [487] 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398 0.883790398
## [493] 0.883790398 0.883790398 0.910092808 0.910092808 0.910092808 0.937699114
## [499] 0.966355040 0.966355040 0.966355040 0.966355040 0.995929592 0.995929592
## [505] 0.995929592 0.995929592 1.026020500 1.057173383 1.057173383 1.057173383
## [511] 1.057173383 1.057173383 1.057173383 1.057173383 1.057173383 1.091579070
## [517] 1.091579070 1.127017733 1.163673474 1.163673474 1.163673474 1.163673474
## [523] 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474 1.163673474
## [529] 1.163673474 1.205390594 1.205390594 1.205390594 1.248789143 1.248789143
## [535] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [541] 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143 1.248789143
## [547] 1.248789143 1.248789143 1.302049650 1.302049650 1.302049650 1.361380792
## [553] 1.423171756 1.423171756 1.423171756 1.423171756 1.423171756 1.492959677
## [559] 1.492959677 1.492959677 1.492959677 1.565134027 1.641566599 1.641566599
## [565] 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599 1.641566599
## [571] 1.641566599 1.641566599 1.760594965 1.760594965 1.760594965 1.760594965
## [577] 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965 1.760594965
## [583] 1.760594965 1.995646102 1.995646102 1.995646102 1.995646102 1.995646102
## [589] 1.995646102 1.995646102 1.995646102 0.000000000 0.000000000 0.000000000
## [595] 0.000000000 0.007077094 0.007077094 0.014216278 0.014216278 0.014216278
## [601] 0.014216278 0.014216278 0.021477375 0.021477375 0.021477375 0.028820277
## [607] 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193 0.036190193
## [613] 0.036190193 0.036190193 0.043683254 0.043683254 0.051238198 0.058848646
## [619] 0.058848646 0.066482062 0.066482062 0.066482062 0.066482062 0.066482062
## [625] 0.066482062 0.066482062 0.066482062 0.066482062 0.074283328 0.082161595
## [631] 0.082161595 0.082161595 0.082161595 0.082161595 0.090312041 0.090312041
## [637] 0.098605212 0.107058768 0.107058768 0.107058768 0.107058768 0.107058768
## [643] 0.115642909 0.115642909 0.124279460 0.124279460 0.124279460 0.124279460
## [649] 0.132971217 0.132971217 0.132971217 0.141833373 0.141833373 0.141833373
## [655] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373
## [661] 0.141833373 0.141833373 0.141833373 0.141833373 0.141833373 0.151197061
## [667] 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061 0.151197061
## [673] 0.151197061 0.151197061 0.151197061 0.160686642 0.160686642 0.160686642
## [679] 0.160686642 0.170396775 0.170396775 0.170396775 0.170396775 0.180213133
## [685] 0.180213133 0.190108901 0.190108901 0.200148605 0.210270329 0.210270329
## [691] 0.210270329 0.210270329 0.210270329 0.220576216 0.220576216 0.220576216
## [697] 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174 0.230965174
## [703] 0.230965174 0.230965174 0.230965174 0.241881805 0.241881805 0.241881805
## [709] 0.241881805 0.241881805 0.241881805 0.252853192 0.252853192 0.252853192
## [715] 0.252853192 0.252853192 0.264161947 0.275576224 0.275576224 0.287173209
## [721] 0.298818473 0.310613892 0.310613892 0.310613892 0.322632099 0.322632099
## [727] 0.322632099 0.322632099 0.322632099 0.334704748 0.334704748 0.334704748
## [733] 0.334704748 0.334704748 0.334704748 0.334704748 0.347024049 0.359487200
## [739] 0.372137832 0.372137832 0.372137832 0.372137832 0.385220199 0.385220199
## [745] 0.398660095 0.398660095 0.398660095 0.398660095 0.398660095 0.412484044
## [751] 0.412484044 0.426505935 0.426505935 0.426505935 0.426505935 0.426505935
## [757] 0.426505935 0.426505935 0.441039991 0.456128884 0.456128884 0.456128884
## [763] 0.471685153 0.471685153 0.471685153 0.471685153 0.487657692 0.487657692
## [769] 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692 0.487657692
## [775] 0.487657692 0.504636398 0.504636398 0.504636398 0.504636398 0.522274168
## [781] 0.540398037 0.540398037 0.558903304 0.577733813 0.597092442 0.597092442
## [787] 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442 0.597092442
## [793] 0.617805328 0.617805328 0.617805328 0.617805328 0.617805328 0.617805328
## [799] 0.617805328 0.617805328 0.617805328 0.640376113 0.640376113 0.640376113
## [805] 0.640376113 0.640376113 0.664790137 0.689971460 0.715491227 0.715491227
## [811] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [817] 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227 0.715491227
## [823] 0.715491227 0.715491227 0.746629486 0.746629486 0.746629486 0.746629486
## [829] 0.746629486 0.746629486 0.780446351 0.780446351 0.780446351 0.815837029
## [835] 0.851972559 0.851972559 0.851972559 0.851972559 0.892199662 0.934636454
## [841] 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454 0.934636454
## [847] 0.934636454 0.934636454 0.934636454 0.934636454 0.986212450 0.986212450
## [853] 0.986212450 0.986212450 1.043419749 1.103396433 1.103396433 1.103396433
## [859] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [865] 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433 1.103396433
## [871] 1.103396433 1.103396433 1.103396433 1.103396433 1.238566987 1.390361656
## [877] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
## [883] 1.390361656 1.390361656 1.390361656 1.390361656 1.390361656
res5<-residuals(Cox.fit5,type="deviance")
hist(res5)
<<<<<<< HEAD
boxplot(res5~Cdata$Death_Censoring)
Opgave 3
boxplot(res5~Cdata$Death_Censoring)
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z p
## factor(Treatment)1 0.067070 1.069371 0.117868 0.569 0.5693
## factor(Sex)1 0.525289 1.690948 0.126974 4.137 3.52e-05
## factor(Ascites)1 0.441588 1.555174 0.179236 2.464 0.0138
## factor(Ascites)2 0.974063 2.648684 0.181719 5.360 8.31e-08
## Age 0.049151 1.050378 0.006764 7.266 3.70e-13
## Prothrombin -0.012923 0.987161 0.002917 -4.431 9.39e-06
##
## Likelihood ratio test=129.5 on 6 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.067070 1.069371 0.117868 0.569 0.5693
## factor(Sex)1 0.525289 1.690948 0.126974 4.137 3.52e-05 ***
## factor(Ascites)1 0.441588 1.555174 0.179236 2.464 0.0138 *
## factor(Ascites)2 0.974063 2.648684 0.181719 5.360 8.31e-08 ***
## Age 0.049151 1.050378 0.006764 7.266 3.70e-13 ***
## Prothrombin -0.012923 0.987161 0.002917 -4.431 9.39e-06 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.0694 0.9351 0.8488 1.3473
## factor(Sex)1 1.6909 0.5914 1.3184 2.1688
## factor(Ascites)1 1.5552 0.6430 1.0945 2.2098
## factor(Ascites)2 2.6487 0.3775 1.8550 3.7819
## Age 1.0504 0.9520 1.0365 1.0644
## Prothrombin 0.9872 1.0130 0.9815 0.9928
##
## Concordance= 0.696 (se = 0.017 )
## Likelihood ratio test= 129.5 on 6 df, p=<2e-16
## Wald test = 132.3 on 6 df, p=<2e-16
## Score (logrank) test = 142.4 on 6 df, p=<2e-16
<<<<<<< HEAD
temp=cox.zph(Cox.fit)
par(mfrow=c(1,2))
plot(temp)
Opgave 4
# pull out different variables to see the effect, leveled varables are not of intresed since all information have be extracted?
# Looking at the model without Age
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit._age <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Prothrombin, data = Cdata)
res_age=residuals(Cox.fit._age)
scatter.smooth(Cdata$Age, res_age)
<<<<<<< HEAD
#Looking at the model without Prothrombin
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit._Prothrombin <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age, data = Cdata)
res_Prothrombin=residuals(Cox.fit._Prothrombin)
scatter.smooth(Cdata$Prothrombin, res_Prothrombin)
<<<<<<< HEAD
Opgave 5
Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
#Interaction Treatment and Sex
Cox.fit_TS <- coxph(Sur1~ factor(Treatment)*factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit_TS
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) * factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z p
## factor(Treatment)1 0.220814 1.247091 0.192829 1.145 0.252157
## factor(Sex)1 0.653601 1.922450 0.181024 3.611 0.000306
## factor(Ascites)1 0.443838 1.558678 0.179358 2.475 0.013339
## factor(Ascites)2 0.967216 2.630611 0.182275 5.306 1.12e-07
## Age 0.049307 1.050543 0.006805 7.245 4.32e-13
## Prothrombin -0.012752 0.987329 0.002927 -4.357 1.32e-05
## factor(Treatment)1:factor(Sex)1 -0.246543 0.781497 0.243803 -1.011 0.311901
##
## Likelihood ratio test=130.6 on 7 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TS)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) * factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.220814 1.247091 0.192829 1.145 0.252157
## factor(Sex)1 0.653601 1.922450 0.181024 3.611 0.000306
## factor(Ascites)1 0.443838 1.558678 0.179358 2.475 0.013339
## factor(Ascites)2 0.967216 2.630611 0.182275 5.306 1.12e-07
## Age 0.049307 1.050543 0.006805 7.245 4.32e-13
## Prothrombin -0.012752 0.987329 0.002927 -4.357 1.32e-05
## factor(Treatment)1:factor(Sex)1 -0.246543 0.781497 0.243803 -1.011 0.311901
##
## factor(Treatment)1
## factor(Sex)1 ***
## factor(Ascites)1 *
## factor(Ascites)2 ***
## Age ***
## Prothrombin ***
## factor(Treatment)1:factor(Sex)1
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.2471 0.8019 0.8546 1.820
## factor(Sex)1 1.9225 0.5202 1.3482 2.741
## factor(Ascites)1 1.5587 0.6416 1.0967 2.215
## factor(Ascites)2 2.6306 0.3801 1.8404 3.760
## Age 1.0505 0.9519 1.0366 1.065
## Prothrombin 0.9873 1.0128 0.9817 0.993
## factor(Treatment)1:factor(Sex)1 0.7815 1.2796 0.4846 1.260
##
## Concordance= 0.695 (se = 0.017 )
## Likelihood ratio test= 130.6 on 7 df, p=<2e-16
## Wald test = 132.1 on 7 df, p=<2e-16
## Score (logrank) test = 142.7 on 7 df, p=<2e-16
#Interaction Treatment and Ascites
Cox.fit_TA <- coxph(Sur1~ factor(Sex)+factor(Treatment)*factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit_TA
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Treatment) * factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z
## factor(Sex)1 0.503528 1.654548 0.126989 3.965
## factor(Treatment)1 0.291419 1.338326 0.139323 2.092
## factor(Ascites)1 0.786840 2.196445 0.234027 3.362
## factor(Ascites)2 1.455465 4.286476 0.243010 5.989
## Age 0.046595 1.047697 0.006718 6.936
## Prothrombin -0.013059 0.987025 0.002968 -4.400
## factor(Treatment)1:factor(Ascites)1 -0.729933 0.481941 0.360030 -2.027
## factor(Treatment)1:factor(Ascites)2 -0.902858 0.405409 0.345164 -2.616
## p
## factor(Sex)1 7.34e-05
## factor(Treatment)1 0.036467
## factor(Ascites)1 0.000773
## factor(Ascites)2 2.11e-09
## Age 4.03e-12
## Prothrombin 1.08e-05
## factor(Treatment)1:factor(Ascites)1 0.042619
## factor(Treatment)1:factor(Ascites)2 0.008904
##
## Likelihood ratio test=139.1 on 8 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TA)
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Treatment) * factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z
## factor(Sex)1 0.503528 1.654548 0.126989 3.965
## factor(Treatment)1 0.291419 1.338326 0.139323 2.092
## factor(Ascites)1 0.786840 2.196445 0.234027 3.362
## factor(Ascites)2 1.455465 4.286476 0.243010 5.989
## Age 0.046595 1.047697 0.006718 6.936
## Prothrombin -0.013059 0.987025 0.002968 -4.400
## factor(Treatment)1:factor(Ascites)1 -0.729933 0.481941 0.360030 -2.027
## factor(Treatment)1:factor(Ascites)2 -0.902858 0.405409 0.345164 -2.616
## Pr(>|z|)
## factor(Sex)1 7.34e-05 ***
## factor(Treatment)1 0.036467 *
## factor(Ascites)1 0.000773 ***
## factor(Ascites)2 2.11e-09 ***
## Age 4.03e-12 ***
## Prothrombin 1.08e-05 ***
## factor(Treatment)1:factor(Ascites)1 0.042619 *
## factor(Treatment)1:factor(Ascites)2 0.008904 **
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1 1.6545 0.6044 1.2900 2.1221
## factor(Treatment)1 1.3383 0.7472 1.0185 1.7585
## factor(Ascites)1 2.1964 0.4553 1.3884 3.4748
## factor(Ascites)2 4.2865 0.2333 2.6623 6.9016
## Age 1.0477 0.9545 1.0340 1.0616
## Prothrombin 0.9870 1.0131 0.9813 0.9928
## factor(Treatment)1:factor(Ascites)1 0.4819 2.0749 0.2380 0.9760
## factor(Treatment)1:factor(Ascites)2 0.4054 2.4666 0.2061 0.7974
##
## Concordance= 0.701 (se = 0.017 )
## Likelihood ratio test= 139.1 on 8 df, p=<2e-16
## Wald test = 148.8 on 8 df, p=<2e-16
## Score (logrank) test = 167.7 on 8 df, p=<2e-16
#standard model
Cox.fit <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z p
## factor(Treatment)1 0.067070 1.069371 0.117868 0.569 0.5693
## factor(Sex)1 0.525289 1.690948 0.126974 4.137 3.52e-05
## factor(Ascites)1 0.441588 1.555174 0.179236 2.464 0.0138
## factor(Ascites)2 0.974063 2.648684 0.181719 5.360 8.31e-08
## Age 0.049151 1.050378 0.006764 7.266 3.70e-13
## Prothrombin -0.012923 0.987161 0.002917 -4.431 9.39e-06
##
## Likelihood ratio test=129.5 on 6 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.067070 1.069371 0.117868 0.569 0.5693
## factor(Sex)1 0.525289 1.690948 0.126974 4.137 3.52e-05 ***
## factor(Ascites)1 0.441588 1.555174 0.179236 2.464 0.0138 *
## factor(Ascites)2 0.974063 2.648684 0.181719 5.360 8.31e-08 ***
## Age 0.049151 1.050378 0.006764 7.266 3.70e-13 ***
## Prothrombin -0.012923 0.987161 0.002917 -4.431 9.39e-06 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.0694 0.9351 0.8488 1.3473
## factor(Sex)1 1.6909 0.5914 1.3184 2.1688
## factor(Ascites)1 1.5552 0.6430 1.0945 2.2098
## factor(Ascites)2 2.6487 0.3775 1.8550 3.7819
## Age 1.0504 0.9520 1.0365 1.0644
## Prothrombin 0.9872 1.0130 0.9815 0.9928
##
## Concordance= 0.696 (se = 0.017 )
## Likelihood ratio test= 129.5 on 6 df, p=<2e-16
## Wald test = 132.3 on 6 df, p=<2e-16
## Score (logrank) test = 142.4 on 6 df, p=<2e-16
<<<<<<< HEAD
# note here that Ascites contains 3 variables, so the p-values can't be interpreted directly so we compair Likelihood ratio tests against the standard model
=======
# note here that Ascites contains 3 variables, so the p-values can't be interpreted directly so we compair Likelihood ratio tests against the standard model
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
#we use a anova co get the chi squared value for the 2 models, the difference in likelihood ratio test is 139.1-129.5=9.4 and the DF is 8-6=2
anova(Cox.fit_TA,Cox.fit) # we get at p-value of 0.008 so it is significant.
## Analysis of Deviance Table
## Cox model: response is Sur1
## Model 1: ~ factor(Sex) + factor(Treatment) * factor(Ascites) + Age + Prothrombin
## Model 2: ~ factor(Treatment) + factor(Sex) + factor(Ascites) + Age + Prothrombin
## loglik Chisq Df P(>|Chi|)
## 1 -1519.4
## 2 -1524.2 9.558 2 0.008404 **
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
#Interaction Treatment and Age
Cox.fit_TAA <- coxph(Sur1~ factor(Sex)+factor(Ascites)+ factor(Treatment)*Age+ Prothrombin, data = Cdata)
Cox.fit_TAA
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + factor(Treatment) *
## Age + Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z p
## factor(Sex)1 0.529397 1.697909 0.126685 4.179 2.93e-05
## factor(Ascites)1 0.462584 1.588172 0.180018 2.570 0.0102
## factor(Ascites)2 0.992521 2.698027 0.182583 5.436 5.45e-08
## factor(Treatment)1 -0.964723 0.381089 0.805313 -1.198 0.2309
## Age 0.040708 1.041548 0.009283 4.385 1.16e-05
## Prothrombin -0.013056 0.987029 0.002919 -4.473 7.71e-06
## factor(Treatment)1:Age 0.016710 1.016851 0.012905 1.295 0.1954
##
## Likelihood ratio test=131.2 on 7 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TAA)
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + factor(Treatment) *
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Sex)1 0.529397 1.697909 0.126685 4.179 2.93e-05 ***
## factor(Ascites)1 0.462584 1.588172 0.180018 2.570 0.0102 *
## factor(Ascites)2 0.992521 2.698027 0.182583 5.436 5.45e-08 ***
## factor(Treatment)1 -0.964723 0.381089 0.805313 -1.198 0.2309
## Age 0.040708 1.041548 0.009283 4.385 1.16e-05 ***
## Prothrombin -0.013056 0.987029 0.002919 -4.473 7.71e-06 ***
## factor(Treatment)1:Age 0.016710 1.016851 0.012905 1.295 0.1954
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1 1.6979 0.5890 1.32458 2.1765
## factor(Ascites)1 1.5882 0.6297 1.11600 2.2601
## factor(Ascites)2 2.6980 0.3706 1.88639 3.8589
## factor(Treatment)1 0.3811 2.6241 0.07862 1.8472
## Age 1.0415 0.9601 1.02277 1.0607
## Prothrombin 0.9870 1.0131 0.98140 0.9927
## factor(Treatment)1:Age 1.0169 0.9834 0.99145 1.0429
##
## Concordance= 0.695 (se = 0.017 )
## Likelihood ratio test= 131.2 on 7 df, p=<2e-16
## Wald test = 134.3 on 7 df, p=<2e-16
## Score (logrank) test = 144.7 on 7 df, p=<2e-16
#Interaction Treatment and Prothrombin
Cox.fit_TP <- coxph(Sur1~ factor(Sex)+factor(Ascites)+ Age+factor(Treatment)*Prothrombin, data = Cdata)
Cox.fit_TP
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + Age +
## factor(Treatment) * Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z p
## factor(Sex)1 0.518523 1.679545 0.127462 4.068 4.74e-05
## factor(Ascites)1 0.441740 1.555411 0.179070 2.467 0.013631
## factor(Ascites)2 0.978996 2.661782 0.181549 5.392 6.95e-08
## Age 0.048780 1.049989 0.006773 7.202 5.93e-13
## factor(Treatment)1 -0.157295 0.854452 0.377810 -0.416 0.677166
## Prothrombin -0.014688 0.985419 0.004080 -3.600 0.000318
## factor(Treatment)1:Prothrombin 0.003450 1.003456 0.005522 0.625 0.532106
##
## Likelihood ratio test=129.9 on 7 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit_TP)
## Call:
## coxph(formula = Sur1 ~ factor(Sex) + factor(Ascites) + Age +
## factor(Treatment) * Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Sex)1 0.518523 1.679545 0.127462 4.068 4.74e-05
## factor(Ascites)1 0.441740 1.555411 0.179070 2.467 0.013631
## factor(Ascites)2 0.978996 2.661782 0.181549 5.392 6.95e-08
## Age 0.048780 1.049989 0.006773 7.202 5.93e-13
## factor(Treatment)1 -0.157295 0.854452 0.377810 -0.416 0.677166
## Prothrombin -0.014688 0.985419 0.004080 -3.600 0.000318
## factor(Treatment)1:Prothrombin 0.003450 1.003456 0.005522 0.625 0.532106
##
## factor(Sex)1 ***
## factor(Ascites)1 *
## factor(Ascites)2 ***
## Age ***
## factor(Treatment)1
## Prothrombin ***
## factor(Treatment)1:Prothrombin
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Sex)1 1.6795 0.5954 1.3083 2.1562
## factor(Ascites)1 1.5554 0.6429 1.0950 2.2094
## factor(Ascites)2 2.6618 0.3757 1.8648 3.7993
## Age 1.0500 0.9524 1.0361 1.0640
## factor(Treatment)1 0.8545 1.1703 0.4075 1.7918
## Prothrombin 0.9854 1.0148 0.9776 0.9933
## factor(Treatment)1:Prothrombin 1.0035 0.9966 0.9927 1.0144
##
## Concordance= 0.695 (se = 0.017 )
## Likelihood ratio test= 129.9 on 7 df, p=<2e-16
## Wald test = 132.4 on 7 df, p=<2e-16
## Score (logrank) test = 142.5 on 7 df, p=<2e-16
<<<<<<< HEAD
Opgave 6
=======Sur1 <- Surv(Cdata$Time, Cdata$Death_Censoring)
Cox.fit <- coxph(Sur1~ factor(Treatment)+ factor(Sex)+ factor(Ascites)+ Age+ Prothrombin, data = Cdata)
Cox.fit
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## coef exp(coef) se(coef) z p
## factor(Treatment)1 0.067070 1.069371 0.117868 0.569 0.5693
## factor(Sex)1 0.525289 1.690948 0.126974 4.137 3.52e-05
## factor(Ascites)1 0.441588 1.555174 0.179236 2.464 0.0138
## factor(Ascites)2 0.974063 2.648684 0.181719 5.360 8.31e-08
## Age 0.049151 1.050378 0.006764 7.266 3.70e-13
## Prothrombin -0.012923 0.987161 0.002917 -4.431 9.39e-06
##
## Likelihood ratio test=129.5 on 6 df, p=< 2.2e-16
## n= 488, number of events= 292
summary(Cox.fit)
## Call:
## coxph(formula = Sur1 ~ factor(Treatment) + factor(Sex) + factor(Ascites) +
## Age + Prothrombin, data = Cdata)
##
## n= 488, number of events= 292
##
## coef exp(coef) se(coef) z Pr(>|z|)
## factor(Treatment)1 0.067070 1.069371 0.117868 0.569 0.5693
## factor(Sex)1 0.525289 1.690948 0.126974 4.137 3.52e-05 ***
## factor(Ascites)1 0.441588 1.555174 0.179236 2.464 0.0138 *
## factor(Ascites)2 0.974063 2.648684 0.181719 5.360 8.31e-08 ***
## Age 0.049151 1.050378 0.006764 7.266 3.70e-13 ***
## Prothrombin -0.012923 0.987161 0.002917 -4.431 9.39e-06 ***
## ---
<<<<<<< HEAD
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
=======
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>>>>>>> e0cec29a27aa2413ceee235b9faad81b4735fba4
##
## exp(coef) exp(-coef) lower .95 upper .95
## factor(Treatment)1 1.0694 0.9351 0.8488 1.3473
## factor(Sex)1 1.6909 0.5914 1.3184 2.1688
## factor(Ascites)1 1.5552 0.6430 1.0945 2.2098
## factor(Ascites)2 2.6487 0.3775 1.8550 3.7819
## Age 1.0504 0.9520 1.0365 1.0644
## Prothrombin 0.9872 1.0130 0.9815 0.9928
##
## Concordance= 0.696 (se = 0.017 )
## Likelihood ratio test= 129.5 on 6 df, p=<2e-16
## Wald test = 132.3 on 6 df, p=<2e-16
## Score (logrank) test = 142.4 on 6 df, p=<2e-16
new.age=(65-57)*365
new.surf<-survfit(Cox.fit,newdata = data.frame(Treatment=0,Sex=1,Age=57,Ascites=1,Prothrombin=85))
autoplot(new.surf,xlab="Days",ylab="Overall survival probability")+
geom_vline(xintercept = new.age)
<<<<<<< HEAD